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Foreword 


This  document  discusses  the  application  of  the  kinematic  wave  routing  method  in  the  Flood 
Hydrograph  Package  (HEC-1)  for  analyzing  urban  runoff  processes.  The  material  is  presented  in  two 
chapters  and  a  comprehensive  example.  Chapter  1  presents  introductory  material  necessary  for  an 
understanding  of  the  theory,  assumptions,  equations  and  numerical  methods  incorporated  into  HEC-1 
for  kinematic  wave  fiood  routing.  The  physical  processes  of  the  urban  runoff  and  streamflow  routing 
are  discussed  briefly  and  related  to  the  kinematic  wave  capabiiities  in  HEC-1 . 

Chapter  2  explains  methods  of  applying  kinematic  wave  routing  techniques  using  HEC-1 .  Data 
requirements  along  with  specific  methods  of  applying  kinematic  wave  routing  techniques  to  runoff 
probiems  in  urban  hydroiogy  are  discussed.  The  chapters  have  been  prepared  so  the  user  can  either 
read  one  or  both.  A  user  interested  oniy  in  the  theory,  or  only  in  application  procedures,  can  read  the 
appropriate  chapter. 

The  basic  approach  for  the  subdivision  of  the  basin  into  the  various  elements  listed  here  was 
deveioped  and  coded  through  discussions  with  John  Peters,  Art  Pabst,  and  Paul  Ely  of  the  Hydologic 
Engineering  Center  (HEC).  The  original  development  and  programming  of  the  kinematic  wave  routing 
routines  in  computer  program  HEC-1  was  done  under  contract  by  Resource  Analysis,  Incorporated, 
under  the  direction  of  Dr.  Brendan  Hariey. 

Subsequently,  the  computational  scheme  was  modified  by  David  Goldman  so  that  the  solution  of 
the  kinematic  wave  equations  is  more  stable  and  provides  the  user  with  information  on  the  accuracy 
of  the  solution. 

Chapter  1  was  originally  written  by  Robert  MacArthur,  and  Chapter  2  was  originaiiy  written  by 
Johnaness  J.  DeVries.  Paul  Detjens  modified  Chapters  1  and  2  to  reflect  the  changes  made  to  the 
program.  He  also  created  a  new  example.  Appendix  A,  which  demonstrates  the  use  of  the  kinematic 
wave  option.  Troy  Nicolini  was  instrumental  in  organizing  this  version  of  TD-10. 
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Chapter  1 

An  Introduction  To  Kinematic 
Fiow  Approximations 


1.1  Introduction 

Rapid  urbanization  of  river  basins  in  and  around  metropolitan  areas  has  forced  land  and  water 
resources  planners  and  hydrologists  to  develop  a  variety  of  methods  for  analyzing  problems  in  urban 
hydrology.  Problems  involving  both  design  and  management  decisions  are  often  so  complex  as  to 
require  application  of  mathematical  models.  The  mathematical  models  which  have  been  most 
commonly  used  rely  on  basic  unit-graph  techniques  to  model  the  application  and  distribution  of 
precipitation  as  rainfail  or  snowmelt,  compute  rainfall  and  snowmelt  losses  and  excesses,  and 
determine  subbasin  outflow  hydrographs.  Although  these  models  which  are  based  on  the 
development  of  a  representative  unit  hydrograph  are  frequently  applied  and  have  been  used 
successfully,  it  is  difficult  to  associate  physical  properties  of  the  basin  to  be  modeled  to  the 
parameters  necessary  to  develop  a  unit  hydrograph.  It  is  even  more  difficult  to  define  some  of  the 
parameters  such  as  the  Clark  and  R  or  Snyder  C^  and  C  for  basins  which  have  no  recorded  data. 
Because  it  is  very  important  to  develop  the  best  representation  of  the  actual  urban  runoff  situation 
when  analyzing  urban  storm  water  runoff  problems,  it  is  desirable  to  relate  runoff  processes  directly  to 
measurable  geographic  features  of  the  basin.  It  is  sdso  desirable  for  the  modeling  technique  to  be 
able  to  reproduce  nonlinear  runoff  characteristics  rather  than  being  limited  to  linear  responses  such  as 
those  developed  by  unit  graph  techniques. 

The  kinematic  wave  method  of  routing  overland  and  river  flows  was  chosen  as  an  additional 
routing  option  for  use  in  the  HEC-1  program  for  several  reasons.  First,  although  simple  in  form, 
kinematic  wave  theory  offers  the  benefits  of  nonlinear  response  without  needing  an  unduly 
complicated  or  costly  solution  procedure.  Second,  for  the  purposes  of  modeling  unsteady  overland 
flow,  any  model  will  require  considerable  parameter  adjustment  to  account  for  the  complexities  of  the 
basin  and  the  specific  flows  which  occur  within  the  basin.  The  kinematic  wave  method  relates  basin 
and  flow  characteristics  directly  to  the  two  routing  parameters,  a  and  m.  The  parameters  a  and  m  are 
directly  related  to  the  shape  of  the  channel,  the  boundary  roughness  and  the  slope  of  the  channel  or 
overland  flow  surface.  There  also  have  been  several  previous  studies  which  have  developed  sets  of 
appropriate  values  for  these  parameters  for  a  large  range  of  flow  and  boundary  conditions.  Third, 
Numerical  techniques  used  to  simulate  overland  and  river  flows  can  only  approximate  the  actual 
response  of  real  systems  because  of  the  complex  nature  of  natural  drainage  basins  and  because 
simplifications  must  be  made  to  the  mathematics  to  make  the  model  efficient  and  economical  to 
execute.  The  kinematic  wave  approximation  has  been  proven  to  be  an  accurate  and  efficient  method 
of  simulating  storm  water  runoff  from  small  basins  for  both  overland  flow  and  stream  channel  routing 
[Overton  and  Meadows,  1976]. 


1.2  Modeling  Unsteady  Flow  Using  the  Kinematic  Wave  Approach 

Kinematics  is  defined  as  the  study  of  motion  exclusive  of  the  influences  of  mass  and  force,  in 
contrast  with  dynamics,  in  which  these  influences  are  included.  Flood  waves  can  be  identified  as 
either  of  two  separate  kinds  of  wave  phenomena:  the  dynamic  wave  and  the  kinematic  wave. 
Although  both  of  these  kinds  of  waves  are  initially  present,  certain  characteristics  of  a  watershed  can 
make  kinematic  waves  the  dominant  characteristic  of  a  flood  event. 


1 


When  inertial  and  pressure  forces  are  important,  ‘dynamic  waves"  govern  the  movement  of  long 
waves  in  shallow  water,  like  a  large  flood  wave  in  a  wide  river  [Stoker,  1957].  When  the  inertial  and 
pressure  forces  are  not  important  to  the  movement  of  the  wave,  ‘kinematic  waves"  govern  the  flow. 

For  this  latter  flow  condition,  the  weight  component  (the  force  in  the  direction  of  the  channel  axis  due 
to  the  weight  of  the  fluid  flowing  downhill  in  response  to  the  action  of  gravity)  is  approximately 
balanced  by  the  resistive  forces  due  to  channel  bed  friction  (in  most  cases  this  is  represented  by 
Manning’s  equation).  Flows  of  this  nature  (kinematic  waves)  will  not  be  accelerating  appreciably  and 
the  flow  will  remain  approximately  uniform  along  the  channel.  No  visible  surface  wave  will  be 
noticeable  and  the  passage  of  the  flood  wave,  as  depicted  in  Figure  1 ,  will  be  seen  by  an  observer  on 
the  bank  as  an  apparently  uniform  rise  and  fall  in  the  water  surface  elevation  over  a  relatively  long 
period  of  time  with  respect  to  the  size  of  the  specific  subbasin  being  analyzed.  Therefore,  kinematic 
flows  are  often  classified  as  uniform  unsteady  flows. 


Groduolly  voried,  unsteody  flow,  streom 
lines  ond  water  surfocc  profiles  ore  not 
porollel. 
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Uniform,  unsteady  flow;  woter  lurfoces  ond 
bed  ore  porollel  to  eoch  other  ond  to  the 
energy  grode  line . 


Figure  1 

Visualization  of  Dynamic  and  Kinematic  Waves 

Dynamic  waves  [Lighthill/Whitham,  1 955]  normally  have  much  higher  velocities  and  attenuate 
more  quickly  than  kinematic  waves,  fiood  flows  are  generally  dominated  by  kinematic  waves.  Even 
though  any  surface  disturbance  will  send  a  "signal"  downstream  at  the  speed  of  small  gravity  waves, 
this  "signal"  will  be  too  weak  to  be  detected  at  any  considerable  distance  downstream.  Therefore,  a 
main  flow  change  or  "signal"  is  carried  as  a  kinematic  wave  at  much  slower  velocities  and  the  speed  of 
flood  waves  may  be  approximated  by  the  speed  of  kinematic  waves.  In  this  context,  a  kinematic  wave 
represents  the  characteristic  changes  in  discharge,  velocity  and  water  surface  elevation  with  time  as 
observed  at  any  one  location  on  an  overland  flow  plane  or  along  a  stream  channel. 
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The  speed  of  small  gravity  waves  occurring  in  shallow  open  channels  is  often  referred  to  as  wave 
celerity,  c  (see  page  vi),  where  the  variable  D  is  the  channel  hydraulic  depth  [Chow,  1959].  The  ratio 
of  the  flow  velocity,  V=Q/A,  to  the  celerity  is  called  the  Froude  number,  F.  Therefore,  the  Froude 
number  also  represents  the  ratio  of  inertial  forces  to  gravity  forces. 

Flows  with  Froude  numbers  greater  than  one  are  classified  as  supercritical  flows  [Chow,  1959], 
and  surface  waves  are  unable  to  move  in  the  upstream  direction  because  the  flow  velocity  is  greater 
than  the  wave  celerity,  c.  Flows  with  Froude  numbers  greater  than  two  tend  to  be  unstable,  which 
may  affect  the  accuracy  and  applicability  of  steady  flow  assumptions  to  these  flows.  However,  the 
characteristics  of  kinematic  waves  [Lighthill/Whitham,  1955]  dominate  over  those  of  dynamic  waves  for 
flows  with  Froude  numbers  that  are  less  than  or  equal  to  two  (e.g.,  F  <  2).  In  fact,  they  found  that  for 
F  <  1,  dynamic  waves  decay  exponentially  with  respect  to  a  time  constant  they  chose  to  define  as 
V/[gS  (1  -  F/2)],  where,  S,  is  the  channel  or  surface  runoff  slope.  Therefore,  one  may  conclude  that 
kinematic  waves  will  ultimately  dominate  the  flow  characteristics  occurring  for  overland  flows  and  small 
watershed  channel  flows  when  the  flow  Froude  number  is  less  than  two. 

This  brief  introduction  to  kinematic  wave  approximations  leads  to  the  following  discussion  and 
definition  of  uniform  unsteady  flows. 


1.3  Review  of  the  Basic  Types  of  Flow 

As  a  brief  review,  several  basic  types  of  open  channel  flow  that  are  most  commonly  experienced 
will  be  discussed  to  remind  the  reader  of  some  of  the  important  differences  and  characteristics  of 
flows  found  in  practice  [Thomas,  1975].  Figure  2  presents  several  possible  types  of  open  channel 
flow. 


UNIFORM 


GRADUALLY  VARIED  ^  ^  ^ 

RAPIDLY  VARIED 


Steady  Flow 


GRADUALLY  VARIED  RAPIDLY  VARIED 


E.L«  •  Energy  line 
W.S.  •  Water  surface 


Unsteady  Flow 


Figure  2 

Possible  Types  of  Open  Channel  Flow  [Thomas,  1975] 
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1.3.1  Steady  Versus  Unsteady  Flow 

If  the  change  in  velcxiity  with  respect  to  time  at  a  given  location  is  zero,  the  flow  is  c^led  steady 
flow.  Otherwise,  the  flow  is  classified  as  unsteady  flow.  Therefore,  unsteady  flows  require  the 
consideration  of  time  as  an  additional  variable. 


1,3.2  Uniform  Venus  Varied  Fiow 


If  the  change  in  channel  flow  velocity,  V,  with  respect  to  distance  along  the  channel,  x,  is  zero 
(i.e.,  dV/dx  =  0)  for  a  given  period  of  time,  the  flow  is  called  uniform.  Otherwise,  the  flow  is 
nonuniform  and  the  relationship  between  kinetic  energy  and  potential  energy  will  be  changing  along 
the  channel.  When  the  flow  is  uniform,  the  water  surface  will  be  parallel  to  the  channel  bottom.  If  the 
flow  is  not  uniform,  the  slope  of  the  water  surface  will  be  slightly  different  than  the  slope  of  the 
channel  bottom. 

If  the  rate  of  change  of  the  water  surface  slope  is  not  visible  to  the  eye,  the  flow  may  be 
considered  as  gradually  varied.  Rapidly  varied  flows  demonstrate  apparent  and  rather  large  water 
surface  slope  changes  (such  as  at  a  hydraulic  jump).  Rapidly  varied  flow  requires  special  treatment 
and  will  not  be  considered  herein. 


1 .4  Development  of  Governing  Equations 

Although  the  basic  differential  equations  capable  of  describing  one  dimensional  gradually  vari^ 
unsteady  flow  were  originally  developed  a  century  ago,  they  have  only  been  recently  applied  (within 
the  last  thirty  to  forty  years)  to  general  hydrologic  engineering  problems  because  it  was  not  possible 
to  solve  these  equations  efficiently  without  a  high  speed  computer. 

The  mechanics  of  unsteady  OF>en  channel  flow  may  be  expressed  mathematically  in  terms  of  the 
equations  developed  in  1870  by  St.  Venant.  Equations  1  and  2  are  partial  differential  equations  which 
may  be  derived  from  the  basic  principles  of  conservation  of  mass  and  momentum.  Various  derivations 
of  these  equations  have  been  presented  [Chow  1959],  [Henderson  1966],  [Strelkoff  1969],  [Fread 
1976]. 


1.4.1  St.  Venant  Equations 


Continuity 


dy 

dt 


(i-f) 


(channels  of  unit  width) 


(1) 


Momentum 


du 

dt 


+  u 


du 

dx 


=  g{S^  -  S^) 


iu-v) 
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where: 


g  =  acceleration  of  gravity  (ft/sec^) 

y  =  mean  depth  (ft) 

q  =  discharge  per  unit  width  of  channel  (cfs/ft) 

X  =  distance  measured  in  downstream  flow  direction  (ft) 

t  =  time  (seconds) 

u  =  x-component  of  mean  velocity  (ft/sec) 

i  =  rainfall  intensity  (cfs/ft^) 

f  =  infiltration  rate  (cfs/ft^) 

Sq  =  average  bottom  slope  (ft/ft) 


(2) 


Sf  =  friction  slope  defined  by  the  Manning  equation 
qL  =  total  lateral  inflow  per  unit  length  of  channel  (cfs/ft^) 

V  =  the  x-component  of  velocity  for  lateral  inflow.  (This  is  assumed  to  be  negligible  to 
the  total  momentum  balance  for  channel  routing  and  is  therefore  zero). 


The  four  terms  in  Equation  1  and  the  five  terms  in  Equation  2,  are  known  successively  as: 
Terms  in  the  Continuity  Equation  (Equation  1) 

•  rate  of  rise  term 

•  storage  term 

•  lateral  inflow  per  unit  length 

•  intensity  of  excess  rainfall 

Terms  in  the  Momentum  Equation  (Equation  2) 

•  acceleration  term 

•  velocity  head  term 

•  depth  taper  term 

•  bed  slope  minus  friction  term 

•  lateral  inflow  term 


Prior  to  presenting  the  details  of  the  kinematic  flow  method  that  have  been  incorporated  into  the 
HEC-1  program,  it  is  important  to  discuss  the  basic  assumptions  and  requirements  associated  with 
the  equations  used  for  gradually  varied  unsteady  flows.  If  these  basic  assumptions  are  not  valid  for 
the  intended  flow  conditions,  then  afternate  methods  of  simulating  the  flow  should  be  sought.  (For 
further  definitions  of  terms,  refer  to  the  Notation  section  on  pages  viii  -  x.) 


1.4.2  General  Assumptions 

In  the  development  of  the  general  unsteady  flow  equations  it  is  assumed  that  the  flow  is  one 
dimensional  in  the  sense  that  flow  characteristics  such  as  depth  and  velocity  are  considered  to  vary 
only  in  the  longitudinal  x-direction  of  the  channel.  Additional  basic  assumptions  necessary  for  the 
validity  of  the  equations  include: 

(1)  the  velocity  is  constant  and  the  water  surface  is  horizontal  across  any  section 
perpendicular  to  the  longitudinal  flow  axis; 

(2)  all  flows  are  gradually  varied  with  hydrostatic  pressure  prevailing  at  all  points  in  the  flow 
such  that  all  vertical  accelerations  within  the  water  column  can  be  neglected; 

(3)  the  longitudinal  axis  of  the  flow  channel  can  be  approximated  by  a  straight  line,  therefore, 
no  lateral  secondary  circulations  occur; 

(4)  the  slope  of  the  channel  bottom  is  small  (less  than  1:10); 

(5)  the  channel  boundaries  may  be  treated  as  fixed  noneroding  and  nonaggrading; 

(6)  resistance  to  flow  may  be  described  by  empirical  resistance  equations  such  as  the 
Manning  or  Chezy  equations; 

(7)  the  flow  is  incompressible  and  homogeneous  in  density. 
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The  assumptions  found  in  items  1  through  7  have  been  shown  [Strelkoff,  1 969]  to  be  valid  and 
applicable  for  most  open  channel  flows  occurring  in  natural  rivers  and  streams.  The  validity  of  the 
restrictions  presented  in  assumption  5,  are  not  easily  evaluated.  Research  is  currently  under  way  to 
estimate  the  overall  affects  of  this  assumption.  It  has  also  been  found  that  overland  flows  such  as 
those  associated  with  stormwater  runoff  can  be  described  by  the  kinematic  wave  form  of  Equations  1 

and  2  wlhoul  violation  of  ttie  seven  previous  assumptions  [Ughtliill/Wtiitham,  1955]  [Liggett/Woolhiser, 
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may  require  special  attention  and  changes  to  methods  currently  used  in  HEC-1 . 


1.5  Discussion  of  the  General  Concepts  and  Structure  of  the  Kinematic  Wave 
Flood  Routing  Technique 

A  general  conceptual  description  of  the  methods  of  modeling  stormwater  runoff  using  kinematic 
wave  simplifications  will  be  presented  here  prior  to  discussing  the  specific  equations  and  algorithms 
that  are  used.  This  will  help  the  reader  understand  how  the  model  ’visualizes"  the  actual  physical 
characteristics  of  the  basin  and  how  it  responds  to  rainfall.  Specific  mathematical  relationships  and 
numerical  solution  techniques  will  be  presented  in  a  following  section. 

Analyses  of  surface  runoff  problems  are  complicated  not  only  by  the  nature  of  specific  storm 
events  which  occur,  but  perhaps  to  a  greater  extent,  by  the  nature  and  complexity  of  the  watershed  or 
urban  area  being  analyzed.  Description  of  the  local  physical  characteristics,  geometry,  and  response 
of  the  system  could  become  a  monumental  task  if  one  were  to  include  every  minute  detail.  The 
concept  incorporated  into  the  HEC-1  program  involves  simulating  the  natural  complexities  of  a  basin 
with  a  number  of  simple  elements,  such  as  overland  flow  planes,  stream  segments,  and  lengths  of 
representative  storm  drain  or  sewer  pipe  as  shown  in  Figures  3  and  4.  If  the  previously  listed 
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Figure  3 

Elements  Used  in  Kinematic  Wave  Caiculations 
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gradually  variad  unstaady  flow  assumptions  ara  not  violatad,  than  combinations  of  thasa  basic 
alamants  hava  provan  to  ba  quita  raprasantativa  of  tha  actual  bahavior  of  most  systams.  Maasurad 
rasponsas  [Wooding,  1966]  from  thraa  natural  drainaga  araas  wara  comparad  with  calculatad  rasults 
obtainad  from  a  modal  which  usad  a  simpla  kinamatic  flow  routing  procadura  similar  to  that  found  in 
HEC-1.  Conclusions  wara  that  although  tha  gaomatry  of  natural  catchmants  is  far  mora  complicatad 
than  that  of  the  simple  model,  the  agreement  between  computed  and  actual  discharge  hydrographs  is 
quite  good. 

To  simulate  the  response  of  a  complex  watershed  to  precipitation  from  storm  events  a 
mathematical  model  made  up  of  combinations  of  simple  geometric  components  is  constructed. 
Successful  application  of  this  approach  begins  with  the  description  of  unsteady,  uniform  flow  over  an 
idealized  planar  overland  flow  element  of  unit  breadth  (called  an  overland  flow  strip)  for  a  number  of 
given  boundary  conditions.  A  boundary  condition  represents  known  or  assumed  flow  conditions  that 
are  specified  by  the  user.  Then,  similar  relationships  for  routing  channel  flows  resulting  from  runoff 
from  overland  flow  elements  are  developed.  Once  these  relationships  have  been  developed, 
combinations  of  simple  elements  can  be  made  to  describe  basin  and  subbasin  responses  to  storm 
events.  Figure  4  summarizes  the  relationship  between  the  three  different  types  of  elements. 


Figure  4 

Relationships  Between  Flow  Elements 


Notice  that  overland  flows  are  handled  separately  from  the  channel  flow,  because  overland  runoff 
demonstrates  specific  shallow  flow  properties  which  lead  to  a  form  of  the  equations  of  motion  which 
differ  from  the  form  required  for  collector  channel  and  main  channel  flows.  These  flows  are  calculated 
individually  and  then  combined  in  such  a  way  as  to  preserve  continuity  and  accuracy.  This 
computational  method  has  been  shown  to  be  quite  efficient  [Harley,  et  al.,  1972]. 
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Figures  3  and  4  present  a  schematic  representation  of  this  approach.  The  governing  equations 
used  for  this  combined  overland  flow  and  river  channel  routing  procedure  are  derived  from  the 
general  St.  Venant  equations  for  unsteady  gradually  varied  open  channel  flow. 

For  the  overland  flow  portion  of  the  model,  shallow  surface  water  runoff  assumptions  are  applied 
to  Equations  1  and  2,  resulting  in  a  simple  kinematic  wave  form  of  the  equations.  Similar 
simplifications  are  found  to  be  valid  and  useful  for  the  channel  routing  portions  as  well.  Attempting 
the  rigorous  mathematical  representation  of  these  complex  phenomena  (runoff  and  channel  routing 
through  the  natural  topography)  would  require  exceedingly  small  spatial  and  temporal  detail  and  result 
in  a  very  large  system  of  simultaneous  equations.  The  kinematic  wave  form  of  the  St.  Venant 
equations  provides  a  simplified  description  of  the  physical  system  in  terms  of  surfaces  and  channels 
with  homogeneous  properties.  The  important  concept  of  the  overland  portion  of  the  model  is  that 
water  is  distributed  over  a  wide  area  and  at  very  shallow  average  depths  until  it  reaches  a  well  defined 
collector  channel.  In  urban  areas,  two  general  types  of  overland  flow  surface  are  usually  present: 
pervious  and  impervious.  The  mechanics  of  flow  over  both  kinds  of  surfaces  are  similar;  however,  the 
slopes,  flow  lengths,  roughnesses,  and  loss  rates  will  differ.  Roof  tops,  parking  lots,  and  paved 
surfaces  such  as  streets,  are  described  as  impervious  areas.  Lawns,  fields,  parks,  etc.,  are  pervious 
areas.  The  percentages  of  the  total  subbasin  area  that  are  impervious  and  pervious  are  stipulated  by 
the  user.  The  model  develops  the  runoff  flows  from  the  rainfall  intensity  and  loss  rates  specified  for 
the  pervious  and  impervious  areas  in  the  basin.  After  the  overland  runoff  is  routed  down  the  length  of 
the  overland  flow  strip  it  is  then  distributed  uniformly  along  the  collector  system  which  represents 
rivulets,  channels,  gutters,  sewers,  and  storm  drains  such  as  shown  in  Figures  3  and  5.  Once  the 
runoff  flows  enter  the  collector  system  (Figure  5)  they  move  through  it  picking  up  additional  (uniformly 
distributed)  lateral  inflow  from  adjacent  runoff  strips.  These  collector  flows  eventually  reach  a  main 
channel  where  they  are  then  routed  as  open  channel  stream  flows  through  the  system. 


Figure  5 

Typical  Urban  Drainage  Pattern 


The  next  sections  will  review  the  basic  concepts  associated  with  the  kinematic  wave  simplification 
of  the  momentum  equation  and  present  some  criteria  [Lighthill/Whitham,  1 955]  that  relates  to  it’s 
vaiidity  for  various  flow  conditions.  Following  sections  will  then  present  the  details  of  the  three 
different  segments  of  flow  routing  (overland,  collector  and  stream  routing)  and  how  they  can  be 
applied  to  problems  in  urban  hydrology. 
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1.6  The  Kinematic  Wave  Form  of  the  Momentum  Equation  is  a  Simpie 
Stage-Discharge  Reiationship 

Recall  that  kinematic  waves  occur  when  the  dynamic  terms  in  the  momentum  equations  are 
negligible.  This  allows  one  to  assume  that  the  bed  slope  is  approximately  equal  to  the  friction  slope 
(Sq  =  Sj).  Under  these  conditions  and  if  there  is  no  appreciable  backwater  effect,  the  discharge  can 
be  described  as  a  function  of  depth  of  flow  only,  for  all  x  and  t. 

Q  =  ay“  (3) 


Q  is  discharge  in  cfs  and  a  and  m  are  kinematic  wave  routing  parameters  which  are  directly  related  to 
the  basin  and  flow  characteristics.  This  can  be  best  appreciated  by  using  an  approach  of 
normalizing  the  momentum  equation  with  a  steady  uniform  discharge  called  Q„  [Henderson,  1966]. 
Rearrangement  of  the  normalized  form  of  Equation  2  yields: 


=  0, 


1 


idy  ^  u  du  1  du 
V  g  dx  g  dt 


+ 


gy ) 


1 

2 


(4) 


If  the  sum  of  the  terms  to  the  right  of  the  minus  sign  is  much  less  than  one  (i.e.,  pressure,  inertia  and 
local  inflow  are  relatively  small  compared  to  SJ,  then  unsteady  flows  are  nearly  uniform  and  may  be 
approximated  by  a  series  of  normal  flows,  e.g., 

Q  “  On  (5) 


Normal  flows  of  this  nature  can  be  described  by  a  depth-discharge  relationship,  such  as  Equation  3 
[Overton/Meadows,  1976].  This  describes  kinematic  flow  and  provides  a  simple  method  for 
calculating  flows  from  stormwater  runoff. 

For  Froude  numbers  less  than  2  [Lighthill/Whitham,  1955],  the  dynamic  component  decays 
exponentially  and  the  kinematic  wave  ultimately  dominates.  As  was  mentioned  in  the  introduction,  this 
means  that  no  visible  surface  wave  is  observed;  only  the  rise  and  fall  of  the  water  surface  can  be  seen 
Figure  1  (page  2).  A  study  [Liggett/Woolhiser,  1967]  on  the  characteristics  of  a  rising  hydrograph  for 
a  large  variety  of  flow  conditions,  has  shown  that  the  dynamic  component  in  Equation  4  will  be 
dampened  enough  to  be  neglected,  provided  that: 

k  =  ^  )  10  (6) 

yF2 

where:  F  =  ^ 

v/gY 

y  =  mean  depth 

=  bed  slope 

k  =  dimensionless  kinematic  flow  number 

L  =  the  length  of  the  plane 

Practical  evaluation  of  k  is  difficult  because  it  may  be  hard  to  estimate  L,  y  or  F  precisely  for  natural 
flow  conditions  without  collecting  field  measurements.  Results  from  the  study  [Liggett/Woolhiser, 

1967]  allow  Equation  2  to  be  greatly  simplified.  These  results  are  summarized  in  Figure  6  (where  Q, 
and  t*  are  dimensionless  discharge;  Q,  =  (qL  L/Vg)  and  time  t*  =  (tV/L),  respectively).  For  a  value  of 
k  of  1 0,  an  approximate  ten  percent  error  in  the  calculated  discharge  hydrograph  would  result  by 
deleting  the  dynamic  terms  from  the  momentum  equation.  Notice  that  as  k  increases  above  ten, 
however,  that  the  error  in  discharge  decreases  rapidly.  A  "true  kinematic’  solution  results  as  k 
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Figure  6 

The  Rising  Hydrograph  -  Variation  with  k,  for  F=1  [LiggettyWoolhiser,  1967] 

approaches  infinity,  but  for  engineering  purposes  flows  characterized  by  k  >  10  can  be  approximated 
reasonably  well  with  the  kinematic  wave  form  of  the  momentum  equation,  e.g..  Equation  3.  Continuing 
research  will  attempt  to  develop  practical  methods  of  evaluating  k  and,  therefore,  the  applicability  of 
kinematic  wave  simplifications  from  typical  basin  data  and  flow  conditions. 

Derivation  of  the  relationships  used  for  overland  flow,  collectors  and  main  channels  will  be 
presented  next. 


1 .7  Development  of  a  and  m  for  Various  Cross  Section  Shapes 

The  following  sections  will  present  example  derivations  of  as  and  m’s  for  the  different  cross- 
sectional  shapes  considered  by  HEC-1 . 


1.7.1  Overland  Flow  Relationship 

The  kinematic  wave  equation  for  an  overland  flow  segment  on  a  wide  plane  with  very  shallow 
flows  can  be  derived  from  Manning's  equation  and  Equation  3.  Consider  the  flow  down  an  overland 
flow  strip  of  unit  breadth  as  shown  in  Figure  3  (page  6).  The  steady  discharge  from  a  unit  strip  can 
be  described  with  Manning's  relationship. 


C  =  idll  rI  sj  A  P) 


where  for  very  shallow  flows  that  are  at  a  depth  of  y^,  R  (the  hydraulic  radius)  and  A  are  simply 
(y^j  •  1)/1  and  y^  •  1)  respectively.  Substitution  of  these  values  for  R  and  A  into  Equation  7,  yields: 


Q  = 


1.486 

N 


(8) 
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Notice  also  that  Manning’s  'n’  has  been  replaced  with  an  appropriate  N  which  describes  the 
properties  of  the  runoff  surface  being  modeled.  Values  of  N  are  usually  greater  than  Manning  s  'n  ’ . 
Table  1  presents  suggested  ranges  of  values  for  N  for  various  surface  conditions.  These  values  were 
obtained  from  previous  field  and  laboratory  investigations  for  overland  flow. 


Table  1 

Effective  Resistance  Parameters  for  Overland  Flow 


Surface 

N  value 

Source 

Asphalt/Concrete* 

0.05  -  0.15 

a 

Bare  Packed  Soil  Free  of  Stone 

0.10 

c 

Fallow  -  No  Residue 

0.008  -  0.012 

b 

Convential  Tillage  -  No  Residue 

0.06  -  0.12 

b 

Convential  Tillage  -  With  Residue 

0.16  -  0.22 

b 

Chisel  Plow  -  No  Residue 

0.06  -  0.12 

b 

Chisel  Plow  -  With  Residue 

0.10-0.16 

b 

Fall  Disking  -  With  Residue 

0.30  -  0.50 

b 

No  Till  -  No  Residue 

0.04  -  0.10 

b 

No  Till  (20  -  40  percent  residue  cover) 

0.07  -  0.17 

b 

No  Till  (60  - 100  percent  residue  cover) 

0.17  -  0.47 

b 

Sparse  Rangeland  with  Debris: 

0.09  -  0.34 

b 

0  Percent  Cover 

0.05  -  0.25 

b 

20  Percent  Cover 

Sparse  Vegetation 

0.053  -  0.13 

f 

Short  Grass  Prairie 

0.10-0.20 

f 

Poor  Grass  Cover  on  Moderately  Rough 

0.30 

c 

Bare  Surface 

Light  Turf 

0.20 

a 

Average  Grass  Cover 

0.4 

c 

Dense  Turf 

0.17  -  0.80 

a,c,e,f 

Dense  Grass 

0.17  -  0.30 

d 

Bermuda  Grass 

0.30  -  0.48 

d 

Dense  Shrubbery  and  Forest  Litter 

0.4 

a 

Legend;  a)  Harley  (1975),  b)  Engman  (1986),  c)  Hathaway  (1945),  d)  Palmer  (1946),  e)  Ragan  and  Duru  (1972),  f)  Woolhlser 
(1975).  [Hjemfelt,  1986] 

Because  the  discharge  represented  by  Equation  8  is  per  unit  breadth,  one  can  substitute  the 
previously  defined  q  for  Q  (see  Equation  2  for  units,  page  4)  to  obtain  discharge  in  terms  of  flow  per 
unit  breadth. 


g  = 


1 . 486  cj  2 


N 


Vo 


(9) 


Rewriting  Equation  3  in  terms  of  discharges  per  unit  width  where  the  subscripts  "o'  indicate  variables 
associated  with  overland  flows. 


•Ho 


Qo  =  “oYo 


(10) 
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where: 

\  N 
=  5/3 

=  average  slope  of  overland  flow  element 
=  mean  depth  for  overland  flow 

Q!q  =  conveyance  for  particular  runoff  surface,  slope,  and  roughness 


Because  there  are  two  unknowns  in  Equation  10  another  relationship  is  required  for  mathematical 
closure  and  a  complete  solution.  A  form  of  the  continuity  equation  (Equation  1 ,  page  4)  provides  the 
necessary  second  equation  to  complete  the  solution. 


dt  dx 


(i  -  f) 


(11) 


where:  (i  -  f)  =  rate  of  excess  rainfall  (rainfall-infiltration)  in  ft/sec 
q^  =  discharge  per  unit  width  in  cfs/ft 
t  =  time  in  seconds 
y  =  mean  depth  of  overland  flow  in  feet 


Together  Equations  10  and  11  form  the  complete  kinematic  wave  equations  for  overland  flow.  If 
Equation  1 0  is  substituted  into  Equation  1 1  one  obtains: 


dt 


(12) 


which  has  only  one  dependent  variable  so  that  it  can  be  solved  to  give  a  relationship  for  in  terms  of 
X,  t,  and  the  excess  rainfall  intensity  (i  -  f).  Once  y^  is  found,  it  can  be  substituted  back  into  Equation 
10  to  obtain  a  value  for  q^.  This  procedure  provides  the  necessary  information  to  be  able  to 
determine  the  time  dependent  discharge  from  the  overland  flow  elements. 


1.7.2  Collector  and  Main  Channel  Routing  Relationships 


For  the  collector  system  (which  represents  rivulets,  storm  drains,  and  sewer  pipes)  and  stream 
channel  segments,  simple  cross  section  shapes  have  been  used  to  simulate  prototype  channels.  It 
has  been  found  that  appropriate  usage  of  simple  triangular,  trapezoidal,  and  circular  channel  shapes 
can  provide  an  accurate  representation  of  the  response  of  the  prototype. 


Flows  entering  the  collectors  and  the  stream  channels  can  consist  of  both  flows  from  upstream 
sections  and  lateral  inflows  from  adjacent  catchment  surfaces.  These  representative  channels  are 
described  by  their  slope,  length,  cross-sectional  dimensions,  shape,  and  Manning's  'n  ’  value.  The 
standard  Manning’s  'n  ’  is  used  here  because  collector  and  streamflows  behave  more  as  normal  open 
channel  flows.  The  basic  form  of  the  equations  for  kinematic  wave  routing  of  collector  and 
streamflows  are  similar  to  those  developed  for  shallow  overland  flow  (Equations  10  and  11).  The 
kinematic  wave  equations  for  collector  and  streamflow  routing  are: 


dt  dx 


(13) 


(14) 
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where: 


Aj,  =  cross-sectional  area  of  flow  in  ft^ 

Qg  =  discharge  in  cfs 

=  lateral  inflow  per  unit  length  in  cfs/ft  from  overland  flow  strips 
t  =  time  in  seconds 

X  =  distance  along  the  stream  in  feet 

Q!g,mj,  =  kinematic  wave  parameters  for  a  particular  cross-sectional  shape,  slope  and 
roughness 


The  reader  should  note  that  the  subscripted  variables  used  above  are  for  a  typical  collector 
channel  and  are  indicated  as  such  with  the  subscript  c.  Identical  relationships  would  be  used  for 
routing  in  the  main  channel  but  one  may  wish  to  identify  them  separately  with  a  subscript  such  as  m 
to  indicate  main  channel  (refer  to  the  Notation  section). 


1 .7.3  Determination  of  and  for  Coliectors  and  Streams 

Values  of  and  m^  will  be  different  for  each  differently-shaped  cross  section  and  will  vary  with 
effective  Manning’s  'n  ’  and  channel  slope  as  well.  The  basic  channel  shapes  considered  by  the 
HEC-1  program  are  the  trapezoid  and  circle.  Variation  of  the  side  slope  and  bottom  width  for  the 
trapezoidal  section  allows  one  to  develop  rectangular  and  triangular  channel  shapes  as  well.  These 
shapes  are  presented  in  Figure  7. 


WITH 
W=  0 


y  TRIANGULAR 

L 


Figure  7 

Two  Basic  Channei  Shapes  and  Their  Variations  Used  by  the  HEC-1 
Program  for  Kinematic  Wave  Stream  Routing 
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As  was  done  for  the  development  of  the  overland  flow  parameters  and  m^,  one  needs  to  define 
the  proper  values  for  and  m^,  for  the  stream  and  collector  system.  Rather  than  derive  all  the  and 
m  values  for  each  differently  shaped  section,  an  easy-to-follow  development  for  a  simple  triangular 
section  will  be  presented  as  an  example.  The  results  for  the  remaining  shapes  will  merely  be 
presented  because  their  derivations  follow  the  same  procedures. 


1.7.4  Triangular  Sections 

Derivation  of  and  m^  for  a  triangular  channel  shape  begins  with  the  description  of  R^,  and  in 
Equation  7  (page  10)  for  a  triangular  shaped  cross  section.  The  reader  should  refer  to  Figure  7  (page 
13)  and  to  the  definition  of  terms  after  Equations  13  and  14. 

=  AJP^  =  (Area  of  triangular  cross  sectionA/Vetted  Perimeter)  =  Hydraulic  Radius 

where:  P,.  =  2  (\Jl  +  tP-  )  •  =  wetted  perimeter 

A^  =  -i  (2yg  •  z)  •  =  zyc^  =  cross-sectional  area 

z  =  side  slope  ratio 


Now  the  values  of  A^,  and  R^.  may  be  substituted  into  Equation  7. 
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Defining  $  as 


1  486  5^ 

,  and  substituting  the  appropriate  value  for  A^  and  into  Equation  1 5 
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and  simplifying,  gives 
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Equation  8  may  now  be  written  in  the  form  of  Equation  14: 


Qc 


=  “c  Ac 


(19) 


where: 


0.94  S/ 
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1.7.5  Rectangular  Sections 

A  rectangular  shape  is  obtained  by  stipulating  that  z  in  Figure  7b  (page  13)  is  zero.  This 
produces  a  channel  w  feet  wide  with  vertical  walls.  Man-made  channels  and  rectangular  concrete 
drain  sections  may  be  represented  by  this  shape.  Following  similar  procedures  which  produced 
Equations  18  and  19,  two  separate  relationships  for  rectangular  shapes  are  easily  derived.  One  for  a 
very  wide  rectangular  channel  where  w  is  much  greater  than  the  depth  y^,  and  the  other  for  a 
rectangular  channel  that  may  have  comparable  depths  and  breadths,  e.g.,  w  ~  y^. 


For  the  wide  shallow  channel  case: 


and 
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(20) 


For  the  rectangular  channel  where  w  «  y^: 


0.7  2  Sc 
n 
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1.7.6  Trapezoidal  Sections 

The  trapezoidal  section  is  one  of  the  two  basic  sections  considered  by  the  general  HEC-1  code. 
Modifications  of  the  trapezoidal  shape  produced  the  rectangular  and  triangular  alternatives  defined 
above.  When  describing  a  trapezoidal  section  it  is  important  to  define  the  side  slopes  z  and  the 
channel  bottom  widths  w  accurately  (Figure  7,  page  13).  It  is  not  possible  to  derive  a  simple 
relationship  for  and  m^.  from  the  geometric  properties  alone,  so  it  becomes  necessary  to  fit  and 
m^,  to  the  Manning  equation  at  two  or  more  depths  y^  and  use  numericai  techniques  of  fitting  the 


“c  = 


and  = 
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kinematic  equation  to  these  values  to  obtain  values  for  and  for  various  flow  conditions. 
Kinematic  wave  equation:  Qc  “  “c  ^ 

Manning  equation  for  a  trapezoid: 


Oa 


1 

1.486  S/ 
n 


_ 1 _ 

w  +  zy^  (1  +  z^)  ^ 


(22) 


where:  A^'  =  the  area  of  the  effective  cross  section  at  depth 


Values  of  m^,  will  vary  from  4/3  for  a  triangular  section  to  5/3  for  a  wide  rectangular  section. 


1.7.7  Circular  Sections 

Circular  sections  can  be  used  to  model  storm  or  sewer  pipes  in  urban  areas.  The  following 
relationships  for  and  m^  were  derived  [Resource  Analysis,  Inc.,  1975]  for  typical  circular  sections 
such  as  that  shown  in  Figure  7a  (page  13)  which  apply  to  pipe  sections  flowing  less  than  90  percent 
full. 


j. 

0-804  i  (23) 

n  ^ 

1.25 

the  diameter  of  the  circular  section  in  feet 


1.8  Numerical  Solution  of  the  Kinematic  Wave  Equations 

The  HEC-1  program  solves  the  kinematic  wave  equations  using  finite  difference  numerical 
techniques.  The  details  of  these  techniques  which  have  been  perfected  after  years  of  development 
and  testing,  are  available  from  several  sources  [Harley,  1975],  [Resource  Analysis,  Inc.,  1975],  [Bras, 
1973].  As  with  standard  numerical  procedures,  time  is  discretized  in  constant  steps  of  At  and 
distance  in  steps  of  Ax.  The  rainfall  excess  (i  -  f)  is  assumed  constant  within  each  time  step  At,  but 
does  change  from  time  step  to  time  step,  to  simulate  the  variability  occurring  within  a  storm  event. 

Recall  that  Equations  13  and  14  were  the  kinematic  wave  equations  for  collector  and  streamflow 
routing.  If  Equation  14  is  substituted  into  Equation  13,  the  following  relationship  which  has  A^.  as  the 
only  dependent  variable  will  be  obtained 


and 


= 

Dc  = 


dt 


(24) 
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Numerical  solution  of  Equation  24  produces  a  relationship  for  in  terms  of  x,  t,  and  q^.  These 
values  of  may  then  be  substituted  into  Equation  14  to  solve  for  Q^.  This  simple  procedure 
(described  previously  for  overland  flows)  allows  the  calculation  of  Q^.  as  a  function  of  the  segment 
length  and  time  t.  Therefore,  one  can  describe  the  discharge  hydrograph  from  each  of  the 
segments  that  are  of  length  L^,.  If  these  discharges,  Q^,  represented  discharges  from  local  collector 
channels,  then  HEC-1  will  distribute  them  uniformly  as  lateral  inflow  into  the  main  channel  or  stream  as 
shown  in  Figures  3  and  4.  Calculation  of  the  resulting  discharge  hydrograph  from  the  main  channel 
or  stream  is  then  performed  in  an  identical  manner  as  was  just  presented.  The  equations  used  would 
be  identical  to  Equations  23  and  24  except  the  subscript  c  would  be  replaced  with  an  m  everywhere 
throughout  Equations  23  and  24.  This  provides  a  simple  straight  forward  procedure  for  calculating 
first  the  overland  flows,  then  the  flows  through  the  collector  channel  system  and  finally  the  discharges 
in  the  main  stream. 

Detailed  development  of  the  specific  finite  difference  equations,  the  coding  procedures  and 
boundary  requirements  can  be  found  in  the  following  references  [Lighthill  and  Whitham,  1955];  [Bras, 
1973]:  [Resource  Analysis,  Inc.,  1975];  [Hydrologic  Engineering  Center,  1979].  The  coding 
modifications  that  were  made  to  include  the  kinematic  wave  routing  procedure  in  HEC-1  was  done  by 
Resource  Analysis,  Inc.,  under  the  direction  of  Dr.  B.  Harley.  The  current  version  of  HEC-1  uses  an 
improved  algorithm  for  solving  the  KW  equations  [HEC,  1990].  The  algorithms  used  are  based  on 
those  developed  for  the  MITCAT  Catchment  Simulation  Model  [Resource  Analysis,  Inc.,  1975].  The 
following  section  presents  a  brief  review  of  the  numerical  methods  and  procedures  utilized  by  the 
HEC-1  program  to  perform  kinematic  wave  routing. 


1.9  Finite  Difference  Solutions  of  the  Kinematic  Wave  Equations 

The  movement  of  a  flood  wave  down  a  stream  or  along  an  overland  flow  surface  can  be  followed 
by  monitoring  the  times  and  locations  where  specific  water  surface  elevations  and  discharges  occur. 
For  example,  if  lateral  inflow  is  assumed  to  be  zero  for  a  moment,  a  flood  wave  can  be  followed  in  time 
by  noting  the  times  when  a  flow  of  a  given  magnitude,  or  a  stage  of  a  given  magnitude  occurs  at 
successive  downstream  stations  along  the  channel.  Flood  routing  methods,  such  as  the  kinematic 
wave  method  used  in  HEC-1,  depend  upon  certain  numerical  techniques  to  solve  the  governing 
equations  which  describe  the  movement,  stage  and  discharge  characteristics  of  a  flood  wave  as  it 
propagates  downstream.  Discussions  of  several  of  the  different  methods  currently  in  use  for 
hydrologic  engineering  studies  are  available  [Overton/Medows,  1976]  [Mahmood/Yevjevich,  1975]. 

The  numerical  method  currently  employed  by  the  HEC-1  program  to  solve  the  governing  equations  is 
a  finite  difference  method.  A  finite  difference  method  (FDM)  presents  a  "pointwise  approximation"  to 
the  governing  partial  difference  equations.  The  FDM  uses  simple  difference  equations  which  replace 
the  partial  differential  equations  for  an  array  of  stationary  grid  points  located  in  the  space-time  (x-t) 
plane  (Figure  8).  The  intersections  of  the  lines  in  Figure  8  define  the  time  and  space  points  at  which 
the  discharge  and  water  surface  elevation  are  computed.  Lines  parallel  to  the  x-axis  are  called  time 
lines,  and  lines  which  are  parallel  to  the  t-axis  are  space  lines.  The  regular  pattern  formed  by  the 
intersections  of  time  and  space  lines  is  called  the  "computational  network."  Those  points  that  are 
marked  by  solid  dots,  represent  computation  points  (called  nodes).  Solutions  to  the  governing 
equations  via  the  FDM  will  be  computed  at  each  of  these  nodes.  Computations  advance  along  the 
downstream  direction  for  each  time  step  At  until  all  the  flows  and  stages  are  calculated  along  the 
entire  distance  L.  Then  the  computation  is  advanced  ahead  in  time  by  one  At  and  the  computations 
for  discharge  and  water  surface  elevation  are  performed  once  again. 

Also  shown  on  Figure  8  are  the  solution  curves  (the  dashed  lines)  that  represent  solutions 
obtained  from  another  method  called  the  "method  of  characteristics."  The  method  of  characteristics  is 
not  an  available  alternate  solution  technique  in  HEC-1 .  As  can  be  seen  in  Figure  8,  the  solution 
curves  (called  characteristic  curves)  do  not  always  intersect  at  a  node  where  evenly  spaced  time  and 
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space  lines  intersect.  Because  of  this,  additional  interpolation  would  be  required  to  obtain  solutions  at 
evenly  spaced  nodes.  To  avoid  this  interpolation,  the  method  of  characteristics  is  not  used  in  HEC-1. 


UPSTREAM  DISTANCE  ^  DOWNSTRE 

STATION  A  STATION  f 


Figure  8 

Characteristic  Curves  on  a  Fixed  Ax  -  At  Grid 


It  is  mentioned  here  though,  because  the  characteristic  curves  represent  locations  in  the  x-t  plane 
where  specific  flow  properties,  such  as  wave  celerity  c,  remain  constant  for  each  time  step  At.  The 
importance  of  this  is  shown  below  when  the  specific  numerical  solution  methods  (the  standard  form 
and  conservation  form  of  the  governing  equations)  are  explained. 

There  are  detailed  methods  available  for  the  approximation  of  derivates  by  finite  differences 
[Carnahan,  et.  al.,  1969].  By  simply  combining  the  appropriate  finite  difference  approximations  for 
first-  and  second-order  derivates,  complete  partial  differential  equations,  such  as  Equations  14  and  24 
can  be  recast  in  terms  of  finite  differences  instead  of  partial  derivatives.  These  new  equations  are 
approximations  of  the  original  equations  but  are  now  in  a  form  which  can  be  easily  handled 
numerically,  especially  with  the  aid  of  high-speed  computers.  As  an  example  of  the  procedure;  the 
first-order  partial  derivative  aQ/3x  is  approximated  using  a  backward  finite  difference  method: 

do  _  Ag  „  ^i.j  ~  (25) 

dx  Ax  Ax 

The  following  sections  will  discuss  how  similar  finite-difference  approximations  are  applied  to  solve  the 
governing  kinematic  wave  equations  in  space  and  time. 
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The  governing  equations  developed  in  the  previous  sections  consisted  of  a  pair  of  equations  for 
each  of  the  different  kinds  of  flow  elements;  e.g.,  Equations  10  and  12  for  overland  flow  elements, 
Equations  13  and  14  for  coliector  elements,  and  two  additional  equations  identical  to  Equations  13 
and  1 4  (only  with  different  subscripts)  for  the  main  channel  routing  elements.  Rather  than  treat  each 
of  these  three  pairs  of  equations  separately,  the  solution  details  will  be  developed  for  one  set  only 
because  they  are  basically  all  the  same.  HEC-1  handles  each  of  these  three  different  kinds  of  flow 
elements  by  using  the  foiiowing  computational  sequence:  computations  start  with  the  determination 
of  overland  flows  which  are  then  input  as  uniformly  distributed  iateral  inflows  into  the  collector 
elements  which,  in  turn,  modify  the  flows  and  distribute  these  collector  flows  uniformly  and  laterally 
along  the  main  channel.  The  main  channel  routes  the  final  flood  wave  through  the  subbasin.  A 
combination  of  several  subbasins  thus  allows  for  the  complete  description  of  an  entire  basin  during  a 
storm  event. 


1.10  Standard  Form  of  the  Kinematic  Wave  Equations 


Consider  Equation  13  which  is  the  continuity  equation  and  Equation  14  which  relates  flow 
within  a  collector  element  to  the  collector  channel  cross-sectional  area  A^.  It  is  assumed  that  the 
kinematic  wave  coefficients  and  m^  are  constant  for  any  given  system  of  channel  elements. 
Differentiation  of  the  flow  Equation  14  with  respect  to  x  and  substitution  of  this  into  Equation  13  gives: 


dt 


Bla-l 

c 


dx 


=  Qo 


(26) 


Notice  that  Equation  26  is  the  same  as  the  previously  derived  Equation  24.  Referring  to  Figure  9,  the 
area  A^  is  known  at  points  on  the  space-time  grid  for  times  prior  to  the  current  time  (designated  by 
the  index  j)  and  at  points  in  space  prior  to  the  current  location  (indicated  by  the  index  i).  Therefore, 
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Figure  9 

Space-Time  Grid  Used  for  Finite  Difference  Method 
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the  index  (i,j)  corresponds  to  the  current  time  and  space  coordinates.  Future  times  and  space 
locations  that  are  advanced  by  one  At  and  Ax  are  indicated  as  j+1  and  i+1,  respectively.  Similarly, 
one  previous  time  and  space  location  would  correspond  to  a  point  on  the  space-time  grid  indexed  by 
(i-1,  j-1).  In  this  way,  each  point  (or  node)  in  the  space-time  grid  can  be  indicated  by  a  double 
subscript  of  i’s  and  j’s.  This  allows  one  to  rewrite  partial  differential  equations  in  terms  of  finite 
difference  approximations  of  known  quantities  (located  at  previous  times  and  space  points)  and 
unknown  quantities  (current  time  and  space  points).  Therefore,  the  area  at  point  B  in  Figure  9 
would  be  designated  as  at  point  C,  ,  j.^j,  at  point  D,  A^^j  and  so  forth.  With  this 

background,  one  can  now  express  the  governing  equations  in  terms  of  finite  differences  using  the 
previously  defined  indexing  scheme  and  solve  for  the  values  of  and  Q^.  at  the  ’current  point  A"  in 
Figure  9. 


amA^  ^ 


djc 


AA 

Ax 


am 


2  Ax 


(27) 


where  the  differential  of  the  area  A  in  the  x-direction  is  taken  as  the  difference  between  the  values 
known  at  points  C  and  D  in  Figure  9.  Also,  the  area  term  which  is  raised  to  the  (m-1)  power  in 
Equation  27  is  considered  to  be  an  average  area  between  points  C  and  D  in  Figure  9. 


Consider  the  time  derivative  term  in  Equation  26.  It  is  evaluated  between  points  A  and  D 

ot 

(Figure  9),  thus: 

dA  AA  _  ( (28) 
dt  At  \  At  ) 

The  lateral  inflow  term  q  is  handled  as  an  average  lateral  inflow  which  occurs  within  a  time  step  At  and 
is  defined  as  q  here  to  simplify  the  final  form  of  the  equation. 


q  «  q(i.j)  q(i.j-i)  =  q 


(29) 


Combining  Equations  27,  28  and  29  produces  the  complete  finite  difference  form  of  the  original  partial 
differential  Equation  26: 
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2  J 

1  Ax 

=  q 


(30) 


Aj  j  is  the  only  unknown  quantity  in  Equation  30  which  can,  therefore,  be  isolated  and  directly  solved 
for 


(m-l) 


(31) 


Once  Ajj  is  known,  the  corresponding  flow  Qj,  can  be  computed  from  Equation  14  (see  page  12). 
This  provides  a  straight  forward  method  of  computing  the  time  varying  discharges  along  the  channel. 
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1.11  Conservation  Form  of  the  Equation 

The  previous  'standard  form"  of  the  equations  applies  in  most  situations  where  the  average  wave 
celerity,  c,  is  less  than  the  ratio  of  the  computational  space  to  the  time  step,  e.g.,  c  <  (Ax/At).  This 
corresponds  to  a  solution  below  a  diagonal  connecting  points  C  and  A  in  Figure  9  (page  19).  When  c 
is  less  than  Ax/At,  this  procedure  provides  an  accurate  approximation  for  the  kinematic  wave 
characteristics  of  a  flood  wave,  and  points  C,  D  and  A  are  used  to  determine  Aj  j. 

However,  if  c  is  greater  than  Ax/At,  is  it  possible  for  flood  wave  characteristics  to  propagate  more 
rapidly  through  space  and  time  than  the  numerical  approximation  method  can  account  for  them.  This 
would  correspond  to  a  solution  found  above  a  diagonal  connecting  A  and  C  in  Figure  9.  (This 
numerical  stability  criteria  can  be  associated  with  the  familiar  Courant  condition  for  stability  of  explicit 
finite  difference  schemes.)  Because  of  the  potential  wave  instability,  an  alternate  form  of  the 
approximate  finite  difference  equations  is  needed. 

The  "conservation  form'  of  the  governing  equations  is  therefore  applied  when  the  average  celerity 
c  of  the  flood  wave  is  greater  than  the  ratio  Ax/At.  In  this  form,  the  temporal  derivatives  are  evaluated 
between  points  B  and  C  in  Figure  9  rather  than  A  and  D,  while  the  spatial  derivatives  are  evaluated 
between  points  B  and  A  rather  than  C  and  D.  In  this  way,  rapidly  advancing  flood  wave 
characteristics  can  be  handled  more  accurately.  HEC-1  checks  the  stability  of  the  wave  for  each  time 
step.  If  it  is  stable,  the  standard  form  of  the  equations  is  used;  if  not,  the  conservation  form  of  the 
equations  is  used. 

The  conservation  form  of  the  spatial  derivative  of  discharge  will  be  evaluated  between  points  B 
and  A  in  Figure  9  (page  19): 

do  _A0  _  Ou.J)  ~  (32) 

dx  Ax  Ax 

and  the  temporal  derivative  of  area  will  be  evaluated  between  points  B  and  C. 

dA  ^  AA  _  ~  (33) 

dt  At  At 

The  substitution  of  these  new  'conservation  form'  derivatives  into  the  continuity  equation  produces: 


Q(i.j)  +  ^(1-1, j)  ■^(1-1, j-i)  _ 


Ax 


At" 


Solving  for  the  only  unknown,  Qj  =  gives 


Ax 


If  desired,  Aj  j  can  be  determined  from  the  flow  equation. 


(34) 


(35) 


(36) 
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1.12  Accuracy  of  the  Finite  Difference  Solution 

The  accuracy  and  stability  of  the  finite  difference  scheme  depends  on  approximately  maintaining 
the  relationship  c(At)=Ax,  where  c  is  the  average  kinematic  wave  speed  in  an  element.  The  kinematic 
wave  speed  is  a  function  of  flow  depth,  and,  consequently,  varies  during  the  routing  of  the  hydrograph 
through  an  element.  Since  Ax  is  a  fixed  value,  the  finite  difference  scheme  utilizes  a  At  variable 
internally  to  maintain  the  desired  relationship  between  ax.  At  and  c.  However,  HEC-1  performs  all 
other  computations  at  a  fixed  time  interval  specified  by  the  user.  Necessarily,  the  variable  At 
hydrograph  computed  for  a  subbasin  by  the  finite  difference  scheme  is  interpolated  to  the  user 
specified  computation  interval  prior  to  other  HEC-1  computations.  The  resulting  interpolation  error  is 
displayed  in  both  intermediary  and  summary  output.  These  interpolation  error  summaries  are  shown 
in  the  example  problem  in  Appendix  A. 

The  accuracy  of  the  finite  difference  scheme  depends  on  the  selection  of  the  distance  increment. 

Ax.  The  distance  increment  is  initially  chosen  by  the  formula  ax  =  cAt,^  where  c,  in  this  instance,  is  an 
estimated  maximum  wave  speed  depending  on  the  lateral  and  upstream  inflows.  At^^  is  the  time  step 
equal  to  the  minimum  of,  (1)  one  third  the  travel  time  through  the  reach,  the  travel  time  being  the 
element  length  divided  by  the  estimated  wave  speed,  (2)  one  fourth  the  upstream  hydrograph  rise 
time,  and  (3)  the  user  specified  computation  interval.  Finally,  Ax  is  chosen  as  the  minimum  of  the 
computed  Ax  and  IVNDXMIN,  where  L  is  the  length  of  the  channel  reach  and  NDXMIN  is  a  user 
specified  number  of  Ax  values  to  be  used  by  the  finite  difference  scheme  (5  <  NDXMIN  <  50  for 
overland  flow  planes,  and  2  <  NDXMIN  <  50  for  channels).  Once  calculated.  Ax  remains  constant  over 
the  routing  through  the  channel  reach. 

Consequently,  the  accuracy  of  the  finite  difference  solution  depends  on  both  the  selection  of  Ax 
and  the  interpolation  of  the  kinematic  wave  hydrograph  to  the  user  specified  computation  interval. 

The  default  selection  of  the  Ax  value  by  the  program  will  probably  result  in  an  error  of  less  than  five 
percent.  Based  on  the  uncertainty  of  the  parameter  assumptions,  this  is  probably  accurate  enough 
for  engineering  purposes.  The  user  may  wish  to  check  the  accuracy  by  altering  NDXMIN  and 
comparing  the  results  (see  Appendix  A).  More  importantly,  the  user  should  always  check  the  error  in 
interpolating  to  the  user  specified  computation  interval  as  summarized  at  the  end  of  the  HEC-1  output. 
This  interpolation  error  may  be  reduced  by  reducing  the  computation  interval,  NMIN. 

This  concludes  the  introductory  development  of  the  general  kinematic  flow  equations  and  a  brief 
discussion  of  their  numerical  solution  techniques.  Chapter  2  will  present  additional  background 
information  and  the  details  necessary  for  the  effective  application  of  these  methods  to  solve  problems 
in  urban  hydrology. 
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Chapter  2 


Application  of  Kinematic  Wave  Routing 
Techniques  using  HEC-1 


2.1  Introduction 

HEC-1  contains  options  for  using  kinematic  wave  theory  to  compute  subbasin  outflow 
hydrographs  and  to  route  hydrographs  through  a  stream  reach  [HEC,  1990].  These  options  provide 
an  alternative  to  the  unit  hydrograph  method  for  determining  direct  runoff.  They  also  provide  a 
streamflow  routing  technique  which  can  be  used  in  place  of  the  Muskingum,  rnodified  Puls,  and  other 
methods  available  in  HEC-1.  As  with  any  type  of  hydrologic  model,  however,  it  is  imperative  th«  the 
modeler  check  the  performance  of  his  modeling  effort  against  observed  data.  Use  of  the  inodel 
without  a  procedure  for  verifying  its  ability  to  correctly  simulate  the  behavior  of  a  given  basin  is 
strongly  discouraged. 

The  purpose  of  this  chapter  is  to  describe  how  the  kinematic  wave  method  in  HEC-1  may  be 
applied.  A  discussion  of  kinematic  wave  theory  as  it  relates  to  this  program  is  giveri  in  Chapter  1  of 
this  document.  Several  papers  are  recommended  for  general  background  on  the  kinematic  wave 
method  [Lighthill/Whitham,  1955],  [Harley,  et  al.,  1972],  [Woolhiser,  1975]. 

One  of  the  attractive  features  of  the  kinematic  wave  approach  to  rainfall-rurioff  modeling  is  that  the 
various  physical  processes  of  the  movement  of  water  over  the  basin  surface,  with  the  attendant 
infiltration,  flow  into  stream  channels,  and  flow  through  the  channel  network  are  considered. 
Parameters  such  as  roughness,  slope,  catchment  length  and  areas,  and  stream  channel  dimensions 
are  used  to  define  the  processes. 

The  various  features  of  the  irregular  surface  geometry  of  the  basin  are  generally  approximated  by 
either  of  two  types  of  basic  elements:  (1)  an  overland  flow  element,  and  (2)  a  stream  or  channel  flow 
element.  In  the  modeling  process  described  here,  one  or  two  overland  flow  elements  (designated  as 
overland  flow  strips)  are  combined  with  one  or  two  channel  flow  elements  to  represent  a  subbasin. 

An  entire  basin  is  modeled  by  linking  the  various  subbasins  together. 

Because  the  descriptions  of  the  various  elements  comprising  the  model  are  directly  related  to 
physical  parameters,  the  model  can  be  easily  modified  to  reflect  changing  land  uses.  This  makes 
kinematic-wave-type  models  very  useful  for  urban  studies  because  the  effects  of  increasing 
urbanization  can  be  accounted  for  by  changing  the  parameters  describing  the  basin. 


The  various  topics  covered  in  this  chapter  include  basin  modeling  procedures,  a  description  of 
the  elements  used  in  kinematic  wave  calculations,  and  procedures  for  selecting  the  parameters.  An 
example  problem,  shown  in  Appendix  A,  is  presented  to  illustrate  HEC-1  input  and  output  data,  and 
the  effects  of  changes  to  numeric  values  of  the  parameters  are  discussed. 
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2.2  Basin  Modeling 


The  modeling  process  starts  with  a  description  of  the  topologic  structure  of  the  basin:  drainage 
basin  boundaries,  stream  and  drainage  channels,  and  the  logical  relationships  between  the  drainage 
areas  and  the  channels.  The  definition  of  the  drainage  boundary  will  depend  on  the  objective  of  the 
study  being  conducted,  as  well  as  the  topological  character  of  the  basin.  Studies  dealing  with  urban 
hydrolo^  usually  require  delineation  of  subbasins  that  are  smaller  than  two  square  miles  in  extent 
(about  five  square  kilometers).  Studies  dealing  with  the  effects  of  channel  modifications  may  permit 
use  of  large  areas;  however,  as  the  area  is  increased  the  assumptions  required  to  apply  the  kinematic 
wave  method  become  more  tenuous.  A  typical  urban  drainage  system  is  shown  in  Figure  5  (page  8) 
Ram  falls  on  two  general  types  of  surfaces:  (1)  those  that  are  essentially  impervious,  such  as  roofs, 
driveways,  parking  lots  and  other  paved  areas;  and  (2)  pervious  areas,  most  of  which  are  covered  with 
vegetation  and  have  numerous  small  depressions  that  produce  local  storage  of  rainfall.  It  is  assumed 
iri  the  model  that  water  initially  travels  over  these  surfaces  as  sheet  flow;  however,  in  a  relatively  short 
distance  the  water  begins  to  collect  in  small  streams  or  rivulets  and  the  process  of  stream  or  channel 
flow  begins.  For  impervious  areas,  the  distance  to  the  first  channel  (e.g.,  a  gutter)  is  typically  thirty  to 
one  hundred  feet.  For  pervious  surfaces,  the  longest  distance  a  drop  of  water  must  travel  to  reach  a 
channel  is  on  the  order  of  one  hundred  to  several  hundred  feet. 


Water  collected  by  the  street  gutters,  travels  no  more  than  a  few  hundred  feet  until  it  enters  catch 
basins  which  are  connected  to  sewers.  These  sewers  are  typically  1.5  to  2  feet  in  diameter  for  the 
local  drains.  The  local  drains  are  connected  in  turn  to  larger  and  larger  drains  which  feed  the  main 
storm  drain.  In  many  areas  the  main  storm  drains  are  open  channels  or  streams.  In  major  urban 
are^,  the  main  storm  drains  are  often  large  closed-conduit  sections,  usually  designed  to  flow  only 

partially  full.  The  kinematic  wave  routing  approach,  which  assumes  open  channel  flow  is  therefore 
appropriate. 


There  are  certain  weaknesses  inherent  in  the  kinematic  wave  routing  approach  which  should  be 
kept  in  mind  by  the  modeler.  These  include  the  following:  (1)  in  kinematic  wave  routing,  the  theory 
does  not  provide  for  attenuation  of  the  flood  wave.  As  a  consequence,  peak  flows  may  be  over 
estirnated.  (2)  Surcharging  of  storm  drains  frequently  occurs  during  major  storm  events,  but  no 
explicit  provision  for  surcharging  is  provided  by  the  method.  (3)  Also,  ponding  and  local  storage  of 
water  during  m^or  events  is  not  accounted  for.  This  might  include  overbank  storage,  ponding  in  the 
streets,  etc.,  which  occur  during  large  storms  of  primary  interest  and  may  not  occur  to  as  great  a 
degree  in  the  events  available  for  calibration  and  verification  of  the  model.  The  modeler  therefore 
should  analyze  the  program  results  to  see  if  this  is  happening. 


2.3  Elements  Used  in  Kinematic  Wave  Calculations 

The  runoff  process  described  above  is  idealized  in  HEC-1  through  the  use  of  the  following  flow 
elements:  (1)  one  or  two  typical  overland  flow  elements,  (2)  a  typical  collector  channel  element,  and 
(3)  a  mam  channel  element.  These  generally  provide  the  necessary  detail  for  modeling  the  runoff 
process  in  urban  basins.  Schematic  drawings  of  these  elements  are  shown  in  Figure  3  (page  6). 
Figure  4  (page  7)  illustrates  the  relationships  between  the  flow  elements.  The  elements  are  specified 
to  represent  typical  features  of  the  basin,  and  thus  the  parameters  chosen  for  the  individual  elements 
should  be  representative  of  the  entire  subbasin.  The  runoff  simulation  process  is  automatically 
expanded  from  the  typical  elements  to  the  whole  subbasin  by  the  program.  Because  land  use  and 
development  practices  are  usually  very  similar  within  a  selected  hydrologic  unit,  assigning  a  single 
value  to  a  given  parameter  usually  gives  good  results. 
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2.4  Overland  Flow  Elements 


The  basic  overland  flow  element  is  simply  a  sloping  rectangular  plane  surface  upon  which  the  rain 
falls.  It  is  modelled  as  a  strip  of  unit  width  (one  foot  or  one  meter  wide).  Some  of  the  rainfall  is  lost  by 
infiltration;  the  remainder  runs  off  the  lower  edge  of  the  plane  into  a  channel.  Infiltration  losses  may 
vary  with  time  or  be  constant,  and  a  different  loss  rate  can  be  specified  for  each  flow  strip.  The 
selection  of  appropriate  loss  rates  is  discussed  in  the  next  section.  The  fraction  of  the  element  that  is 
impervious  can  also  be  specified. 

The  basic  kinematic-wave-analysis  concept  used  in  HEC-1  allows  the  use  of  either  one  or  two 
overland  flow  surfaces,  each  discharging  into  a  collector  channel.  For  example,  one  element  could 
represent  all  areas  that  are  essentially  impervious,  with  short  lengths  of  flow  (LJ  to  the  point  where 
the  flow  becomes  channel  flow.  This  element  would  represent  driveways,  roofs,  street  surfaces,  etc. 

The  other  overland  flow  element  could  represent  areas  that  are  pervious  and  have  higher 
resistance  to  flow,  such  as  lawns,  fields,  and  wooded  areas.  In  general,  the  catchment  flow  lengths 
and  roughness  coefficients  will  be  much  greater  for  these  areas.  Again,  the  value  of  to  be  used  is 

the  representative  maximum  distance  for  water  to  travel  as  overland  flow  for  this  type  land  surface. 

The  user  of  this  method  should  think  of  the  overland  flow  strips  as  representing  typical  flow 
surfaces  rather  than  actual  planar  surfaces,  except  when  very  small  areas  (such  as  one  city  lot)  are 
being  considered.  It  is  only  at  these  very  small  scales  that  the  mean  surface  slope  and  actual  area 
and  length  come  close  to  fitting  the  basic  theoretical  concept. 

The  following  data  are  needed  as  input  to  HEC-1  to  describe  each  overland  flow  strip: 

-  typical  overland  flow  length 

-  representative  slope 

N  -  roughness  coefficient  (Table  1,  page  11) 

Aoi.Ao2  -  The  percentages  of  the  subbasin  area  which  the  overland  flow  surface  represents 
(two  possible  types  for  each  subbasin) 


Methods  to  determine  these  parameters  are  discussed  in  the  following  sections.  It  is  suggested  that 
the  data  first  be  tabulated  on  a  form  for  the  entire  basin  and  then  entered  into  an  input  data  file  for 
HEC-1.  A  typical  data  tabulation  sheet  is  given  in  Table  2. 


2.5  Loss  Rate  Selection 

Loss  rates  are  a  very  important  component  of  the  kinematic  wave  model,  and  must  be 
specified  for  the  overland  flow  strips.  Unfortunately,  these  loss  rates  are  often  difficult  to  exactly 
determine.  When  the  analysis  focuses  on  an  ungaged  basin,  as  is  often  the  case  in  urban  hydrology, 
the  SCS  curve  number  method  is  the  most  popular.  This  curve  number  technique  [SCS,  1975] 
involves  estimating  a  curve  number  from  soil  hydrologic  groups  identified  on  a  soil  survey  map.  The 
SCS  has  recommended  using  the  curve  numbers  shown  in  Table  3  for  use  in  urban  areas. 

Many  of  these  curve  numbers  for  urban  land  use  are  based  on  percent  impervious  data.  The 
curve  numbers  presented  here  are  actually  weighted  averages  using  a  curve  number  of  98  for 
impervious  regions,  and  a  value  between  39  and  80  for  the  pervious  areas.  For  the  kinematic  wave 
model,  the  pervious  loss  rate  should  be  selected  based  on  pervious  area,  such  as  open  space,  not  on 
the  land  use  under  study.  If  the  pervious  curve  number  is  selected  based  on  land  use,  the  effect  of 
the  percent  impervious  will  be  counted  twice  and  result  in  an  overestimation  of  runoff. 
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Table  2 

Parameters  for  Kinematic  Wave  Example  Problem 


Percent 

Total 

Su  basin 
Area 

L(ft) 

S 

Ftoughness 

Chanttel 

Shape 

Channel 
Size  (ft) 

Z 

Collector 

System 

Area 

(ml®) 

Loss 

Rate 

(In/hr) 

SubbMin  1  (1 .5  ml^ 

■ 

1 .  Overland  Flow  Strip 

1CX) 

500 

0.04 

0.40 

- 

- 

H 

- 

0.50 

2.  Collector  Channel 

- 

1,500 

0.025 

0.10 

Triangle 

- 

1 

0.40 

- 

3.  Main  Channel 

- 

3,500 

0.010 

0.05 

Triangle 

“ 

4 

- 

- 

Subbaain  2  (1 .2  ml^ 

■ 

1 .  Overland  Flow  Strip  1 

20 

50 

0.06 

0.30 

- 

- 

H 

- 

0.02 

2.  Overland  Flow  Strip  2 

80 

180 

0.01 

0.40 

- 

- 

H 

- 

0.20 

3.  Collector  Channel 

- 

2,100 

0.008 

0.020 

Triangle 

- 

1 

0.35 

“ 

4.  Main  Channel 

- 

4,000 

0.003 

0.025 

Trapezoid 

2 

2 

- 

- 

An  example  will  help  clarify  this  point: 

For  both  unit  hydrograph  and  kinematic  wave  models,  determine  the  correct  loss  rate  for  a  subbasin  with  Vi  acre 
residential  lots  on  a  soil  group  of  A.  The  land  use  in  the  subbasin  is  considered  to  be  uniform. 

•  Untt  Hydrograph  Modal;  the  curve  number  corresponding  to  Vi  acre  lots  on  soil  group  A  is  61  (Table  3).  This 
would  be  the  correct  loss  rate  to  use  for  the  subbasin. 

•  Klftematic  Wava  Modal:  This  can  be  modeled  using  two  overland  flow  stripe,  one  representing  impervious  areas 
(CN=90),  the  other  representing  the  pervious  areas,  such  as  open  space  and  lawns.  A  curve  number  of  39  would 
be  appropriate  for  these  pervious  areas  on  soil  group  A.  The  average  1/4  acre  lot  Is  38%  Impervious  (Table  3), 
which  means  38%  of  the  basin  area  would  have  a  curve  number  of  98,  and  the  remaining  62%  would  have  a  curve 
number  of  39.  Thus,  98  and  39  would  be  the  appropriate  loss  rates  to  use  on  the  two  overland  flow  strips. 

Note  that  with  this  data,  a  composite  curve  number  can  be  calculated:  (0.38)  *  98  +  (0.62)  *  39  =  61 

Since  this  composite  curve  number  (used  for  the  unit  hydrograph  model)  includes  the  percent  impervious,  it  should 
not  be  used  to  model  the  pervious  areas  in  the  kinematic  wave  model.  Using  this  curve  number  would  result  in  an 
overestimation  of  the  impervious  runoff. 

Another  concern  in  the  selection  of  ioss  rates  is  the  proper  consideration  of  unconnected 
impervious  areas.  In  developed  areas  of  the  watershed,  most  impervious  areas  are  directly  connected 
to  storm  drain  systems.  However,  some  impervious  areas  may  drain  onto  pervious  surfaces  which 
reduces  the  volume  of  the  resulting  runoff.  For  example,  the  part  of  a  roof  which  drains  into  a 
backyard  or  landscaping  would  be  considered  as  unconnected  impervious  area  because  the  total 
volume  of  the  resulting  runoff  would  probably  not  enter  the  storm  drain  system.  The  effects  of 
unconnected  imperviousness  should  be  considered  when  determining  loss  rates.  Site  inspection  or 
sensitivity  analysis  could  be  used  to  determine  the  extent  of  unconnected  imperviousness.  Section 
2.6  discusses  how  to  correct  for  unconnected  imperviousness. 
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Table  3 

Runoff  Curve  Numbers  for  Selected  Agricultural, 
Suburban,  and  Urban  Land  Use 


Antecendent  Moisutre  Condition  II.  and  1^  =  0.25  [SCS,1975] 


Hydrologic  Soil  Group 

Land  Use  Description 

D 

D 

C 

!D 

Cultivated  Land^ :  without  conservation  treatment 

72 

81 

88 

91 

with  conservation  treatment 

62 

71 

78 

81 

Pasture  or  Range  Land:  poor  condition 

68 

79 

86 

89 

good  condition 

39 

61 

74 

80 

Meadow:  good  condition 

30 

58 

71 

78 

Wood  or  Forest  Land:  thin  stand,  poor  cover,  no  mulch 
good  cover^ 

45 

66 

77 

83 

25 

55 

70 

77 

Open  Spaces,  Lawns,  Parks,  Golf  Courses,  Cemeteries,  etc. 

good  condition:  grass  cover  on  75  percent  or  more  of  the  area 

39 

61 

74 

80 

fair  condition:  grass  cover  on  50  to  75  percent  of  the  area 

49 

69 

79 

84 

Commercial  and  Business  Areas  (85  percent  impervious) 

89 

92 

94 

95 

Industrial  Districts  (72  percent  impervious) 

81 

88 

91 

93 

Residential:^ 

Average  Lot  Size  Average  Percentage  Impervious^ 

acre  or  less  65 

77 

85 

90 

92 

14  acre  38 

61 

75 

83 

87 

14  acre  30 

57 

72 

81 

86 

14  acre  25 

54 

70 

80 

85 

1  acre  20 

51 

68 

79 

84 

Paved  Parking  Lots,  Roofs,  Driveways,  etc.® 

98 

98 

98 

98 

Streets  and  Roads: 

Paved  with  curbs  and  storm  sewers® 

98 

98 

98 

98 

Gravel 

76 

85 

89 

91 

Dirt 

72 

82 

87 

89 

^  For  a  more  detailed  description  of  agricultural  land  use  curve  numbers  refer  to  National  Engineering  Handbook,  Section  4, 
Hydrology,  Chapter  9,  August  1972. 

^  Good  cover  is  protected  from  grazing  and  litter  and  brush  cover  soil. 

®  Curve  numbers  are  computed  assuming  the  runoff  from  the  house  and  driveway  is  directed  towards  the  street  with  a 
minimum  of  roof  water  directed  to  lawns  where  additional  infiltration  could  occur. 

^  The  remaining  pervious  areas  (lawn)  are  considered  to  be  in  good  pasture  condition  for  these  curve  numbers. 

®  In  some  warmer  climates  of  the  country  a  curve  number  of  95  may  be  used. 
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2.6  Correction  for  Unconnected  Impervious 

The  effects  of  unconnected  impen/ious  can  be  accounted  for  by  using  a  composite  curve  number. 
This  curve  number  (CN)  is  a  function  of  the  pervious  area  CN,  the  percent  impervious  area,  and  the 
percent  of  impervious  area  that  is  unconnected.  Figure  10  shows  a  graph  which  can  be  used  to 
determine  the  composite  curve  number  [McCuen,  1989]. 

On  the  right  hand  side  of  the  graph,  find  the  intersection  of  the  total  percent  impervious  and  the 
percent  unconnected.  Move  horizontally  to  the  left  from  this  intersection  until  the  intersection  with  the 
line  representing  the  pervious  (not  composite)  CN.  The  number  directly  below  this  second 
intersection  is  the  desired  composite  CN. 

This  adjustment  is  not  without  limitations.  For  the  composite  CN  to  be  valid,  the  total  percent 
impervious  must  be  less  than  30  percent.  If  the  total  percent  impervious  is  too  large,  the  additional 
runoff  from  the  unconnected  areas  cannot  be  assumed  to  be  uniformly  distributed  on  the  pervious 
area.  Also,  as  the  percent  impervious  increases,  the  remaining  pervious  regions  become  quickly 
saturated  and  thus  less  effective  in  infiltrating  water. 


Composite  CN 


Total  Impervious  Area  (%) 


Figure  1 0 

Graph  for  Determining  Composite  Curve  Numbers 


To  use  this  composite  curve  number  in  a  kinematic  wave  model,  it  must  be  broken  down  into  it’s 
pervious  and  impervious  components.  The  impervious  CN  remains  98  and  the  total  percent 
impervious  is  assumed  to  be  the  same,  so  the  following  equation  can  be  used  to  calculate  the 
adjusted  pervious  curve  number. 


X  = 


CN^-98*f 

iTf 


(37) 


where:  X 


f 


adjusted  pervious  curve  number 
composite  curve  number  (Figure  1 0) 
total  fraction  impervious  0^<1 
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This  adjusted  pervious  curve  number  can  then  be  used  as  the  loss  rate  for  the  pervious  overland  flow 
strip.  A  curve  number  of  98  can  be  used  for  the  impervious  flow  strip. 

2.7  Selection  of  Overland  Flow  Parameters 

The  area  is  the  simplest  quantity  to  specify:  the  area  of  each  element  is  giveti  as  a  percentage 
of  the  total  area  of  the  subbasin.  If  a  single  element  is  used,  one  hundred  percent  is  sprecified.  If  two 
elements  are  used,  the  sum  of  the  percentages  should  be  one  hundred  percent. 

The  slope  is  a  value  which  is  representative  of  the  slope  of  the  path  that  the  water  takes  on  its 
way  to  the  collector  channel.  It  may  differ  from  the  mean  topographic  slope  for  the  catchment,  and  it 
is  usually  strongly  related  to  the  type  of  land  use  or  development.  For  an  urban  setting,  a  single 
slope-value  can  be  used  for  all  areas  of  similar  building  practice,  even  though  the  mean  ground 
slopes  vary  significantly. 

As  discussed  in  Chapter  1 ,  the  kinematic  wave  equations  for  the  wide  flow  planes  of  the 
overland  flow  elements  are  based  on  Manning’s  equation  for  flow  in  a  wide  channel: 

mo  (38) 

qo  =  “o  Yo 


where:  q^ 

y 


flow  per  unit  width 
flow  depth 

kinematic  wave  parameter 
kinematic  wave  parameter 


For  this  situation: 


=  5/3  and  =  (1 .49/N)S''' 


where:  S  =  slope 

N  =  surface  roughness  coefficient  similar  to  Manning’s  'n  ’  for  channel  flow 


Because  the  nature  of  sheet  flow  with  very  small  depths  over  rough  surfaces  differs  markedly  from 
streamflow,  these  resistance  factors  have  much  different  values  than  the  Manning’s  'n  ’  values  used  in 
streamflow’ computations.  Values  of  N  found  to  be  appropriate  are  given  in  Table  4.  Use  of 
Manning’s  'n  ’  instead  of  these  roughness  coefficients  will  result  in  exaggerated  values  for  q^^  on  the 
overland  flow  strips. 

A  critical  parameter  in  the  overland  flow  element  description  is  the  flow  length  L^.  A  proper 
specification  of  is  vital,  since  it  is  the  most  important  parameter  in  determining  the  response 
characteristics  oAhe  overland  flow  elements.  The  overland  flow  length  can  be  looked  on  as  the 
maximum  length  of  the  path  taken  by  a  representative  water  drop  to  reach  a  channel  where  it  first 
moves  as  streamflow.  It  is  thus  the  distance  for  overland  flow  to  reach  a  tributary  or  local  channel, 
such  as  the  street  gutter. 

Fortunately,  in  many  natural  basins  and  urban  catchments,  close  examination  of  the  full  drainage 
system  reveals  that  the  small-scale  drainage  patterns  are  quite  similar  throughout  the  entire  basin. 
The  value  of  L  appropriate  for  such  a  situation  will  not  vary  greatly  over  the  basin.  However,  the 
actual  values  A  which  give  the  correct  runoff  response  for  the  basin  must  be  verified  through 
comparison  of  model  output  with  measured  data. 


31 


Table  4 

Effective  Resistance  Parameters  for  Overland  Flow 


Surface 

N  value 

Source 

Asphalt/Concrete* 

0.05  -  0.15 

a 

Bare  Packed  Soil  Free  of  Stone 

0.10 

c 

Fallow  -  No  Residue 

0.008  -  0.012 

b 

Convential  Tillage  -  No  Residue 

0.06  -  0.12 

b 

Convential  Tillage  -  With  Residue 

0.16  -  0.22 

b 

Chisel  Plow  -  No  Residue 

0.06  -  0.12 

b 

Chisel  Plow  -  With  Residue 

0.10  -  0.16 

b 

Fall  Disking  -  With  Residue 

0.30  -  0.50 

b 

No  Till  -  No  Residue 

0.04-0.10 

b 

No  Till  (20  -  40  percent  residue  cover) 

0.07-0.17 

b 

No  Till  (60-100  percent  residue  cover) 

0.17  -  0.47 

b 

Sparse  Rangeland  with  Debris: 

0.09  -  0.34 

b 

0  Percent  Cover 

0.05  -  0.25 

b 

20  Percent  Cover 

Sparse  Vegetation 

0.053  -  0.13 

f 

Short  Grass  Prairie 

0.10  -  0.20 

f 

Poor  Grass  Cover  on  Moderately  Rough 

0.30 

c 

Bare  Surface 

Light  Turf 

0.20 

a 

Average  Grass  Cover 

0.4 

c 

Dense  Turf 

0.17-0.80 

a,c,e,f 

Dense  Grass 

0.17-0.30 

d 

Bermuda  Grass 

0.30  -  0.48 

d 

Dense  Shrubbery  and  Forest  Litter 

0.4 

a 

^Asphalt/Concrete  n  value  for  open  channel  flow  0.01  -  0.016 

Legend:  ^  Hadey  (1975),  b)  Engman  (1986),  c)  Hathaway  (1945),  d)  PaJmer  (1946),  e)  Ragan  and  Duru  (1972),  1)  Woolhiser 
(1975).  [Hjemfeft,  1986] 


2.8  Collector  Channel 

The  collector  channel  element  is  used  to  model  the  flow  in  its  path  from  the  point  where  it  first 
becomes  channel  flow  to  the  point  where  it  enters  the  main  channel.  The  inflow  to  the  collector 
channel  is  taken  as  a  uniformly  distributed  flow  along  the  entire  length  of  the  channel.  This  correctly 
represents  the  situation  where  overland  flow  runs  directly  into  the  gutter,  and  also  provides  a 
reasonable  approximation  of  the  flow  inputs  into  the  storm  drain  system  from  individual  catch  basins 
and  tributary  collector  pipes  distributed  along  the  collector  channel. 

The  kinematic  wave  equations,  developed  in  Chapter  1,  are  the  continuity  equation  for  unsteady 
channel  flow  with  lateral  inflow,  and  Manning’s  equation. 


The  continuity  equation  is: 


dt 


(39) 


where: 

Ac 

% 

X 


=  the  channel  flow  in  cfs 
=  flow  cross-sectional  area  in  square  feet 
=  lateral  inflow  in  cfs/ft  to  the  channel 

=  distance  along  the  channel  In  feet  and  t  is  the  time  in  seconds 


Manning’s  equation  is  written  in  the  form: 


Qc 


IDfi 

ttoAc 


(40) 


Where  the  kinematic  wave  routing  coefficients  a  and  m  are  a  function  of  the  channel  geometry.  The 
genereil  expression  for  a  is: 


a 


C 


Ky/ST 

n 


(41) 


where: 


K 

Sc 


n 


=  a  constant  that  depends  on  the  channel  geometry 
=  channel  slope 

=  Manning’s  roughness  coefficient 


A  change  in  either  n  or  S^,  will  change  the  value  of  used  in  the  calculations.  It  should  be  noted  that 
the  model  is  more  sensitive  to  changes  in  n  than  to  S^,  because  depends  on  the  first  power  of  n 
while  it  is  proportional  to  the  square  root  of  S^. 

The  value  of  the  exponent  m^  for  trapezoidal  channels  ranges  from  4/3  when  the  trapezoid  has  a 
base  width  of  zero  (triangular  shape)  to  5/3  for  a  very  wide  rectangular  shape.  For  a  channel  with  a 
circular  cross  section,  m^  is  taken  as  1.25.  See  Chapter  1  for  the  derivation  of  these  parameters. 

The  following  data  are  needed  as  input  to  describe  the  collector  channel  system  (refer  to 
Figure  7,  page  13): 

•  The  surface  area  drained  by  a  single  representative  collector  channel  (e.g.,  gutter  plus  storm 
drain),  A,,. 

•  The  collector  channel  length  (total  length  of  gutter  plus  length  of  storm  drain),  L^. 

•  The  channel  shape  (either  a  circular  section  or  some  variant  of  a  trapezoid). 

•  The  pipe  diameter  or  the  trapezoid  bottom  width  and  side  slope,  of  appropriate. 

•  The  channel  slope,  S^. 

•  Manning’s ‘n’. 


33 


2.9  Selection  of  Collector  Channel  Parameters 


The  characteristics  of  the  collector  channel  components  can  be  determined  by  looking  at  a 
drainage  map  for  the  basin  and  selecting  a  typical  collector  system  for  each  subbasin.  This  single 
typical  collector  channel  is  used  to  represent  all  the  collector  channels  in  the  subbasin. 

The  area  associated  with  the  collector  system  can  be  determined  from  the  map.  This  is  an  area 
in  square  miles  (or  square  km)  rather  than  the  percentage  of  the  subbasin  area.  It  does  not  have  to 
be  an  integer  multiple  of  the  subbasin  area. 

The  collector  channel  length  is  taken  as  the  longest  flow  path  from  the  upstream  end  of  the 
collector  system  to  its  outlet  at  the  main  channel.  This  length  should  include  the  distance  the  water 
will  travel  as  gutter  flow. 

The  channel  shape  and  size  will  usually  change  along  the  length  of  the  channel;  however,  a 
single  shape  must  be  chosen  to  represent  the  channel  along  its  entire  length.  This  is  not  as  great  a 
problem  as  might  appear  at  first.  As  shown  in  Figure  11,  a  triangle  with  side  slope  of  to  1:1  matches 
reasonably  well  with  the  area-discharge  relationships  for  circular  conduits  for  a  given  slope  and 
roughness.  HEC-1  assumes  a  triangular  shape  if  the  channel  shape  is  not  specified  in  the  input  data. 
The  selection  of  channel  shape  is  discussed  in  more  detail  in  Section  2.10. 


Figure  1 1 

Discharge  Versus  Flow  Area  for  Various  Cross  Sections 


If  the  representation  is  by  a  circular  or  trapezoidal  shape,  the  chosen  channel  dimensions  should 
represent  the  most  commonly  used  size  of  channel  in  the  system. 

The  channel  slope  can  be  estimated  from  a  topographic  map  by  taking  the  difference  in  elevation 
between  the  upstream  and  downstream  ends  and  dividing  by  the  length.  If  drop  structures  are  used 
in  the  storm  drains,  the  slopes  should  be  adjusted  accordingly. 
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A  Manning’s  'n  ’  which  best  represents  the  roughness  of  the  major  portion  of  the  channel  should 
be  used.  Tables  of  'n  ’  for  various  types  of  channels,  such  as  concrete  pipe,  line  or  unlined  open 
channels,  etc.,  are  available  in  hydraulic  handbooks  and  other  sources. 


2.10  Main  Channel 

The  main  channel  can  carry  both  inflows  from  upstream  subbasins  as  well  as  flows  supplied  by 
the  collector  channels  within  the  subbasin.  The  inflow  from  the  collector  channel  is  taken  to  be 
uniformly  distributed  along  the  length  of  the  main  channel.  This  is  assumed  to  reasonably 
approximate  the  actual  situation  where  the  flow  enters  the  channel  from  various  collectors  at  a  number 
of  discrete  points  at  various  spacings.  In  Equation  39,  the  lateral  inflow,  q^,  is  determined  by  scaling 
up  the  collector  channel  flow  to  match  the  total  subbasin  areas  and  then  dividing  the  flow  by  the  total 
main  channel  length. 


An  example  will  help  clarify  this  point; 

Suppose  that  the  subbasin  area  is  1 .0  square  mile,  while  the  collector  channel  area  is  0.30  square  mile.  T^e  length  of 
the  main  channel  is  2,000  feet  If  the  collector  channel  flow  is  designated  as  Q^,  the  inflow  per  foot  of  mam  channel  is: 


a  -  /  t-O  [cjfs/ft] 

^0.30^2000^ 

The  channel  routing  element  can  also  be  used  independently  for  routing  a  hydrc^raph  through  a 
channel  reach.  If  desired,  the  subbasin  flow  can  be  computed  separately  and  combined  with  a  routed 
flow  at  the  subbasin  outlet.  Any  of  the  routing  methods  available  in  HEC-1  can  be  used  for  channel 
routing  (Muskingum,  modified  Puls,  Tatum,  etc.)  if  desired. 


The  channel  routine  procedure  requires  the  following  data: 


•  Channel  or  stream  length, 

•  Slope, 

•  Manning’s  'n  ’ 

•  Area  of  subbasin,  A^^bbasin 

•  Channel  shape  (trapezoidal  or  circular) 

•  Channel  dimensions  (e.g.,  width  (w),  side  slopes  (z),  or  diameter  (D^)) 

•  The  upstream  hydrograph  to  be  routed  through  the  reach  if  desired. 


2.11  Selection  of  Main  Channel  Parameters 

Most  of  the  channel  data  can  be  obtained  from  physically  measurable  parameters  for  the  channel 
and  subbasin,  with  the  exception  of  Manning’s  'n  ’ .  The  following  procedure  can  be  used  to 
determine  these  data: 

1 .  The  channel  length  can  be  scaled  from  a  drainage  map  of  the  basin. 

2.  The  mean  channel  slope  can  be  obtained  from  field  measurements  or  estimated  using 
topographic  maps. 

3.  Selection  of  Manning’s  *n  ’  should  be  based  on  the  average  channel  conditions. 

4.  The  subbasin  area  can  be  measured  from  topographic  maps. 

5.  The  selection  of  channel  cross  section  is  discussed  in  Section  2.12. 
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6.  The  channel  dimensions  follow  from  the  preceding  item. 

7.  An  upstream  hydrograph  will  not  be  routed  through  the  channel  reach  unless  the  user 
specifically  requests  the  HEC-1  program  to  do  so  (RK  Record,  Field  8) 


2.12  Selection  of  Channel  Cross  Section 

The  channel  cross  section  for  either  the  main  channel  or  the  collector  channel  can  be  defined  as 
one  of  two  simple  shapes  to  permit  modeling  a  variety  of  natural  channels.  The  kinematic  wave  model 
is  not  especially  sensitive  to  channel  cross-sectional  shape  in  the  simulation  of  discharge,  and 
therefore,  it  is  not  necessary  to  use  complex  channel  shapes.  The  shapes  which  can  be  used  in 
HEC-1  are  trapezoidal,  deep  rectangular  and  circular.  Most  main  channels  can  be  best  simulated  by 
a  triangular  shape.  For  a  main  channel,  this  can  be  done  by  specifying  the  base  width  for  the 
trapezoid  as  zero.  For  a  collector  channel,  the  trapezoid  base  width  defaults  to  zero.  When  the 
channel  is  small,  flood  flows  generally  require  overbank  areas  to  carry  the  flows,  and  a  triangular 
shape  usually  represents  this  situation  quite  well.  As  seen  by  a  plot  of  area  versus  flow.  Figure  1 1 ,  a 
triangle  can  be  used  to  approximate  circular  sections  as  well  as  street  gutters. 

In  the  downstream  reaches  of  the  basin,  the  trapezoidal  section  is  usually  the  most  appropriate 
choice  for  open  channels.  The  chosen  shape  should  provide  the  best  fit  to  the  channel  shape  for  the 
given  flood  flow,  however.  For  example,  in  some  situations  a  trapezoid  might  be  the  best  fit  for  low 
flows,  while  a  triangle  might  be  more  appropriate  for  high  flows.  An  illustration  of  this  is  shown  in 
Figure  12. 


b )  Representotion  by  a  Trapezoid 


Mailmum  Flow 
D«pth 


c )  Representotion  by  a  Triongte 


Figure  1 2 

Fitting  Channel  Cross-Sectional  Shape 
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The  circular  section  allows  modeling  of  storm  sewers.  The  flow  behavior  of  the  conduit  is 
simulated  properly  up  to  a  point  where  the  conduit  is  approximately  ninety  percent  full.  The  program 
does  not  handle  the  effects  of  pressure  flow,  so  for  flows  greater  than  about  ninety  percent  of  pipe 
capacity,  HEC-1  assumes  that  the  capacity  of  the  element  increases  as  required  and  has  no  upper 
limit.  In  many  cases  this  approach  adequately  represents  what  is  happening  in  nature,  because  water 
that  does  not  enter  the  storm  drains  flows  over  the  surface  or  along  a  roadway  until  it  finds  another 
location  to  enter  the  drainage  system. 


2.13  Summary  and  Conclusions 

The  material  presented  in  this  chapter  provides  the  user  of  HEC-1  with  background  information  for 
modeling  watersheds  basins  using  kinematic  wave  routing.  The  example  problem  in  Appendix  A  gives 
an  illustration  of  the  use  of  the  method  and  highlights  the  internal  processes  used  in  HEC-1  to 
minimize  errors.  Guidance  on  selection  of  parameters  is  provided,  along  with  a  brief  sensitivity 
analysis  to  show  effects  of  varying  the  parameters. 

It  is  important  for  the  user  of  this  method  to  verify  results  by  comparing  them  to  measured  rainfall- 
runoff  events.  This  comparison  assesses  the  performance  of  the  model  with  the  selected  parameters. 
Without  such  a  check,  the  inexperienced  modeler  should  interpret  the  results  with  a  great  deal  of 
caution.  However,  kinematic  wave  models  have  been  used  successfully  in  urban  hydrology  in  a  large 
number  of  applications,  and  when  the  models  are  properly  formulated,  good  simulations  result. 
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Appendix  A 

An  Example  Application  of  Kinematic  Wave  Methods 


The  small  partially-urban  basin  shown  in  Figure  A-1  is  to  be  modeled  using  kinematic  wave  runoff 
and  routing  options  of  HEC-1 .  The  hydrologic  characteristics  of  this  basin  (obtained  by  previous 
calibration)  are  as  foilows: 

Subbasin  1 :  The  upper  of  the  two  subareas  making  up  the  basin  is  presently  not  urbanized  and 
is  primarily  rolling  pasture  land  with  few  trees.  The  typical  distance  for  flow  to  travel  to  tributary 
stream  channels  is  500  feet.  The  overland  flow  roughness  coefficient,  N,  is  0.4.  The  representative 
ground  slope  is  0.04.  The  amount  of  impervious  area  is  assumed  to  be  negligible,  and  the 
subbasin  area  A^.,  is  1.5  square  miles. 


Figure  A>1 

Basin  for  Example  Problem 


The  collector  or  tributary  channels  have  a  typical  slope  S^.  of  0.025  and  an  ‘n’ value  of  0.10,  with 
a  typical  channel  length,  L^.,  of  1,500  feet.  The  most  representative  section  is  a  triangle.  The  area, 

Ag ,  contributing  to  a  typical  collector  stream  is  0.4  square  mile.  The  main  channel  is  approximately 
triangular  in  cross  section  with  side  slopes,  z,  of  1  in  4.  The  mean  channel  slope  is  0.01  and 
Manning’s  'n  ’  is  0.05.  Its  length,  L^,  is  3500  feet. 

Subbasin  2;  The  lower  subbasin  is  completely  urbanized.  Twenty  percent  of  the  subbasin 
surfaces  are  impervious.  In  this  subbasin  the  impervious  runoff  areas  have  the  following 
characteristics:  =  50  feet,  =  0.06,  N  =  0.15.  The  pervious  areas  can  be  represented  by  the 

following  parameters:  =  130  feet,  =  0  °‘''  ^  =  0.30.  The  subbasin  area,  A^g,  is  1.2  square 

miles.  The  total  basin  area  is  2.7  square  miles.  The  collector  channel  system  involves  three  hundred 
feet  of  gutter  plus  an  additional  1,800  feet  of  pipe  storm  drain  ranging  up  to  four  feet  in  diameter.  A 
triangular  section  is  used  to  represent  the  various  channel  components  (the  program  default  value 
with  one  to  one  side  slopes).  The  average  slope,  Sp,  is  0.008,  and  the  Manning’s  'n’, which  accounts 
for  friction  and  other  channel  head  losses  is  0.020.  The  area.  A,  contributing  to  the  collector  channel 
system  is  0.35  square  mile. 
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The  parameters  describing  this  basin  are  given  in  Table  A-1 .  Also,  a  listing  of  the  program  input 
for  modeling  this  basin  is  provided  in  Figure  A-2  and  hydrographs  from  the  run  output  are  shown  in 
Figure  A-3.  A  listing  of  the  program  output  is  provided  in  Figure  A-4. 

Table  A-1 

Parameters  for  Kinematic  Wave  Example  Problem 


Percent 
Total 
SubbasI 
n  Area 

L(ft) 

S 

Roughness 

Channel 

Shape 

Channel 
Size  (ft) 

Z 

Collector 

System 

Area 

(ml®) 

Loss 

Rate 

(In/hr) 

Subbasin  1  (1 .5  ml^ 

■ 

■ 

1 .  Overiand  Flow  Strip 

too 

500 

0.04 

0.40 

Hi 

■ 

0.50 

2.  Collector  Channel 

- 

1,500 

0.025 

0.10 

- 

1 

0.40 

- 

3.  Main  Channel 

- 

3,500 

0.010 

0.05 

- 

4 

- 

- 

Subbasin  2  (1  ^  ml^ 

■ 

1 .  Overiand  Flow  Strip  1 

20 

50 

0.06 

0.30 

_ 

_ 

__ 

_ 

2.  Overiand  Flow  Strip  2 

80 

180 

0.01 

0.40 

- 

- 

- 

- 

3.  Collector  Channel 

- 

2,100 

0.008 

0.020 

Triangle 

- 

1 

0.35 

mm 

4.  Main  Channel 

- 

4,000 

0.003 

0.025 

Trapezoid 

2 

2 

- 

■■ 

10 

Kinematic  Wave  Routing  In  HEC-1 

10 

Example  Basin  Model  Used  In  Training  Document  10 

10 

10 

Original  Basin  Conditions 

IT 

5 

1MAY79 

1200  100 

10 

3 

KK 

SUBI 

KM 

Subarea  1 

A=1.5  Sq  Ml  L=3500  ft 

KO 

1 

0 

PB 

2,0 

ZR 

=P1  A= 

TEST  B=EXAMPLE  C=PREC1P-1NC  F=06S 

LU 

0.0 

0.5 

0.0 

BA 

1  .5 

UK 

500 

.04 

.4  100 

RK 

1500 

.025 

.1  .4 

RK 

3500 

.01 

.05  2.17  TRAP  0  4 

ZW 

A=EXAMPLE  C=FL0W  F=C0MPUTED 

KK 

SUB2 

KM 

Subaea  2 

A=1.2  Sq  Ml  L=400  ft 

KO 

1 

0 

LU 

0.0 

0.02 

0.0  0.0  0.20  0.0 

BA 

1 .2 

UK 

50 

.06 

.3  20 

UK 

180 

.01 

.4  80 

RK 

2100 

.008 

.02  .35 

RK 

4000 

.003 

.025  0  TRAP  2  2  YES 

ZW 

C=FL0W 

ZZ 

Figure  A-2 

HEC-1  Input  File  Using  Kinematic  Wave  Method 
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FLOOD  HYDROGRAPH  PACKAGE  (HEC-I)  • 
SEPTEMBER  1090  • 

VERSION  4.0  • 

RUN  DATE  04/30/1001  TIME  14:32:40  • 


U.S.  ARMY  CORPS  OF  ENGINEERS 
HYDROLOGIC  ENGINEERING  CENTER 
600  SECOND  STREET 
DAVIS,  CALIFORNIA  05616 
(016)  756-1104 


X  X 

X  X 

X  X 

XXXXXXX 


xxxxxxx 

X 

X 

xxxx 

X 

x 

xxxxxxx 


xxxxx 

X  X 

X 
X 
X 

X  X 

xxxxx 


X 

XX 

X 

xxxxx  X 
X 
X 

XXX 


THIS  PROGRAM  REPLACES  ALL  PREVIOUS  VERSIONS  OF  HEC-1  KNOWN  AS  HEC1  (JAN  73),  HECIGS,  HECIOB,  AND  HEC1KW. 


THE  KFINITIONS  OF  VARIABLES  -RTIHP-  ANO  -RTIOR-  HAVE  CHANGED  FROM  THOSE  USED  WITH  THE  1B73-STVLE  INPUT  fiTnilrTinc 

i:?{LlRS^^S 


LINE 


HEC-1  IW>UT 

10 . 1 . 2 . 3 . 4 . 5....... 8....... 7... _ 8, 


0 . 10 


PAGE  1 


1 

2 

a 

4 

5 

6 

7 

8 
0 

10 

11 

12 

13 

14 

15 

16 


ID 

ID 

10 

10 

IT 

10 

KK 

KM 

KO 

PB 

PI 

LU 

BA 

UK 

RK 

RK 


Kin«Mtlc  Wave  Routing  In  HEC-1 
Exaaplo  Batin  Modal  Utad  in  Training  Doeuaant  10 


Original  Batin  Conditlont 
5  1UAY76  1200  100 

4 


SUBI 

Subarea  1 

A-1.5 

Sq  Mi  L«:3500  ft 

3 

2.0 

6 

2 

12 

15 

25  20 

0.0 

0.5 

0.0 

1.5 

500 

.04 

.4 

100 

1500 

.025 

.1 

.4 

3500 

.01 

.06 

2.17  TRAP 

17 

KX 

SUB2 

16 

KM 

Subaaa  2 

A-1.2  8q 

Mi  LMOO  ft 

10 

KO 

1 

2 

20 

LU 

0.0 

0.02 

0.0 

0.0 

0.20 

21 

BA 

1.2 

22 

UK 

50 

.06 

.3 

20 

23 

UK 

160 

.01 

.4 

80 

24 

RK 

2100 

.008 

.02 

.35 

25 

RK 

4000 

.003 

.025 

0 

TRAP 

26 

ZZ 

FLOOD  HYDROGRAPH  PACKAGE  (HEC-1)  • 

SEPTEMBER  1900  • 

VERSION  4.0  • 

RUN  DATE  04/30/1001  TIME  14:32:40  • 


2  YES 


U.S.  ARMY  CORPS  OF  ENGINEERS 

•  •• 

HYDROLOGIC  ENGINEERING  CENTER 

609  SECOND  STREET 

DAVIS,  CALIFORNIA  05616 

(016)  756-1104 

Kinautlc  Wave  Routing  in  HEC-1 
Exaapla  Bttln  Modal  Utad  in  Training  Oocuaant  10 

Original  Batin  Conditlont 


6  10 


IT 


OUTPUT  CONTROL 

VARIABLES 

IPRNT 

4 

PRINT  CONTROL 

IPLOT 

0 

PLOT  CONTROL 

QSCAL 

0. 

HYDROGRAPH  PLOT  SCALE 

HYDROGRAPH  TIME 

DATA 

NMIN 

5 

MINUTES  IN  COMPUTATION  INTERVAL 

IDATE 

1MAY79 

STARTING  DATE 

ITIME 

1200 

STARTING  TIME 

NO 

100 

NUMBER  OF  HYOROGRAPH  ORDINATES 

NOOATE 

1MAY79 

ENDING  DATE 

NOTIME 

2015 

ENDING  TIME 

ICENT 

10 

CENTURY  MARK 

COMPUTATION  INTERVAL 

.08  HOURS 

TOTAL  TIME  BASE 

8.25  HOURS 

Figure  A-4 

HEC-1  Output  File  for  Example  Problem 
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ENGLISH  UNITS 

DRAINAGE  AREA 
PRECIPITATION  DEPTH 
LENGTH,  ELEVATION 
FLOW 

STORAGE  VOLUME 
SURFACE  AREA 
TEMPERATURE 


SQUARE  MILES 

INCHES 

FEET 

CUBIC  FEET  PER  SECOND 

ACRE -FEET 

ACRES 

DEGREES  FAHRENHEIT 


7  KK  •  SUB1 


Subarea  1  A«1.5  8q  Ml  L«3500  ft 

9  KO  OUTPUT  CONTROL  VARIABLES 

IPRNT  3  PRINT  CONTROL 

IPLOT  2  PLOT  CONTROL 

QSCAL  0.  HYDROGRAPH  PLOT  SCALE 


SUBBASIN  RUNOFF  DATA 

13  BA  SUBBASIN  CHARACTERISTICS 

TAREA  1.50  SUBBASIN  AREA 

PRECIPITATION  DATA 

10  P6  STORM  2.00  BASIN  TOTAL  PRECIPITATION 

11  PI  INCREMENTAL  PRECIPITATION  PATTERN 

6.00  12.00  15.00  25.00  20.00  10.00  6.00  4.00 


UNIFORM  LOSS  RATE 

STRTL 

.00 

INITIAL 

LOSS 

CNSTL 

.50 

UNIFORM 

LOSS  RATE 

RTIMP 

.00 

PERCENT 

IMPERVIOUS  AREA 

14  UK 


KINEMATIC  WAVE 

OVERLAND-FLOW  ELEMENT 
L  500. 

8  .0400 

N  .400 

PA  100.0 

OXMIN  5 


NO.  1 

OVERLAND  FLOW  LENGTH 
SLOPE 

ROUGHNESS  COEFFICIENT 

PERCENT  OF  SUBBASIN 

MINIMUM  NUMBER  OF  OX  INTERVALS 


15  RK 


16  RK 


KINEMATIC  WAVE 
COLLECTOR  CHANNEL 


L 

1500. 

CHANNEL  LENGTH 

S 

.0250 

SLOPE 

N 

.100 

CHANNEL  ROUGHNESS  COEFFICIENT 

CA 

.40 

CONTRIBUTING  AREA 

SHAPE 

TRAP 

CHANNEL  SHAPE 

WO 

.00 

BOTTOM  WIDTH  OR  DIAMETER 

2 

1 .00 

SIDE  SLOPE 

NOXMIN 

2 

MINIMUM  NUMBER  OF  OX  INTERVALS 

MAIN  CHANNEL 

L 

3500. 

CHANNEL  LENGTH 

8 

.0100 

SLOPE 

N 

.050 

CHANNEL  ROUGHNESS  COEFFICIENT 

CA 

1.50 

CONTRIBUTING  AREA 

SHAPE 

TRAP 

CHANNEL  SHAPE 

WD 

.00 

BOTTOM  WIDTH  OH  DIAMETER 

2 

4.00 

SIDE  SLOPE 

NOXMIN 

2 

MINIMUM  NUMBER  OF  OX  INTERVALS 

RUPSTQ 

NO 

ROUTE  UPSTREAM  HYDROGRAPH 

COMPUTED  KINEMATIC  PARAMETERS 
VARIABLE  TIME  STEP 
(OT  SHOWN  IS  A  MINIMUM) 


ELEMENT 

ALPHA 

M 

OT 

OX 

PEAK 

TIME  TO 

VOLUME 

MAXIMUM 

PEAK 

CELERITY 

(MIN) 

(FT) 

(CFS) 

(MIN) 

(IN) 

(FPS) 

PLANE1 

.75 

1 .67 

5,01 

100.00 

2250.61 

36.15 

1.63 

.33 

COLLECTOR 1 

1.18 

1 .33 

1 .21 

500,00 

2181 .44 

40.61 

1.63 

8.02 

MAIN 

1.16 

1 .33 

1.61 

1166.67 

2151 .60 

41 .02 

1.63 

10.82 

CONTINUITY  SUMMARY  (AC-FT)  -  INFLOW-  ,OOOOE+O0  EXCESS-  .1333E+03  OUTFLOW-  .1304E+03  BASIN  STORAGE-  . 36896+00  PERCENT  ERROR-  1.9 

INTERPOLATED  TO  SPECIFIED  COMPUTATION  INTERVAL 


MAIN 

1 . 

16  1 

.33  5,00 

2137.53 

... 

• 

•  • 

•** 

... 

HYDROGRAPH 

AT  STATION  SU61 

TOTAL  RAINFALL  - 

2.00,  TOTAL 

LOSS  « 

.33,  TOTAL  EXCESS  » 

1.67 

PEAK  FLOW 

TIME 

MAXIMUM  AVERAGE 

FLOW 

(CFS) 

(HR) 

e-HR 

24 -HR 

72 -HR 

e.25-HR 

2138. 

.67 

(CFS) 

261 . 

191 . 

191, 

191  . 

(INCHES) 

1 ,620 

1 .630 

1.630 

1.630 

(AC-FT) 

130, 

130. 

130. 

130. 

CUMULATIVE 

AREA  » 

1.50  SO  MI 

Figure  A-4  (continued) 
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STATION 


8uei 


OAHRMN 
11200 
11205 
11210 
11215 
11220 
11225 
11230 
11235 
11240 
11245 
11250 
11255 
11300 
11305 
11310 
11315 
11320 
11325 
11330 
11335 
11340 
11345 
11350 
11355 
11400 
11405 
11410 
11415 
11420 
11425 
11430 
11435 
11440 
11445 
11450 
11455 
11500 
11505 
11510 
11515 
11520 
11525 
11530 
11535 
11540 
11545 
11550 
11655 
11600 
11605 
11610 
11615 
11620 
11626 
11630 
11635 
11640 
11645 
11650 
11655 
11700 
11705 
11710 
11715 
11720 
11725 
11730 
11735 
11740 
11745 
11750 
11755 
11800 
11805 
11810 
11815 
11820 
11825 
11830 
11835 
11840 
1 1845 
11850 
11855 
11  GOO 
11905 
11910 
11015 
11920 
11925 
11930 
11935 
11940 
11045 
11950 
11055 
12000 
12005 
12010 
12015 


0. 

.0 

PER 

10--- 

20 

3.0 

4.  0 

5. 

6. 

7. 


0. 
10. 
11 . 
12, 

13. 

14. 

15. 

16. 

17. 

18. 

19. 

20. 
21 . 
22. 

23. 
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19  KO  OUTPUT  CONTROL  VARIABLES 
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CELERITY 

(MIN) 

(FT) 

(CFS) 

(MIN) 

(IN) 

(FPS) 

PLANE1 

1.22 

1 .67 

.62 
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INTERPOLATED  TO  SPECIFIED  COMPUTATION  INTERVAL 
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RUNOFF  SUNUARY 

FLOW  IN  CUBIC  FEET  PER  SECOND 
TIME  IN  HOURS,  AREA  IN  SQUARE  MILES 


OPERATION 

STATION 

PEAK 

FLOW 

TIME  OF 
PEAK 

AVERAGE 

6<KOUR 

FLOW  FOR  MAXIMUM  PERIOD 
24-HOUR  72-HOUR 

BASIN 

AREA 

MAXIMUM 

STAGE 

TIME  OF 
MAX  STAGE 

HYDROGRAPH 

AT 

SUBl 

2138. 

.67 

261 . 

1B1. 

191 . 

1 .50 

HYDROGRAPH 

AT 

SUB2 

4030. 

.67 

500. 

365. 

365. 

2.70 

ISTAQ  ELEMENT 


DT 
(UIN) 

1,91  2151.60 


SUMMARY  OF  KINEMATIC  WAVE  -  MUSKINGUH-CUNGE  ROUTING 
{FLOW  IS  DIRECT  RUNOFF  WITHOUT  BASE  FLOW) 

INTERPOLATED  TO 

^  COMPUTATION  INTERVAL 

PEAK  TIME  TO  VOLlAiE  DT  PEAK  TIME  TO 

PEAK  PEAK 

(CFS)  {MIN)  (IN)  (MIN)  (CF6)  (MIN) 


VOLUME 

(IN) 

1.63 


6UB1  MANE  1,91  2151.60  41,02  1.63  5.00  2137.53  40.00 

CONTINUITY  SUMMARY  (AC-FT)  •  INFLOW=  .OOOOE+00  EXCESS-  .1333Et03  OUTFLOW-  .1304E*03  BASIN  STORAGE-  .3689E-f00  PERCENT  ERROR-  1.9 
SUB2  MANE  1.47  4086.93  37.47  1.73  5.00  4030.04  40.00  1.73 

CONTINUITY  SUMMARY  (AC-FT)  -  INFLOW-  .1304E+03  EXCESS-  .1210Et03  OUTFLOW-  .24B3Et03  BASIN  STORAGE-  .2237E+00  PERCENT  ERROR-  .7 


***  NORMAL  END  OF  HEC-1  ••• 


Figure  A-4  (continued) 


The  table  entitled  'SUMMARY  OF  KINEMATIC  WAVE  -  MUSKINGUM-CUNGE  ROUTING',  found  at  the  end  of 
the  output  listing,  is  especially  important  for  kinematic  wave  modeling.  The  user  should  always  examine  the 
difference  between  the  peak  discharge,  time  to  peak  and  hydrograph  volume  computed  with  an  internal  At  with 
that  interpolated  to  the  user  specified  At. 

Consider,  for  example,  the  results  summarized  for  subbasin  2  in  Figure  A^.  The  computed  peak  flow  is 
4,086.93  rfs,  the  time  to  peak  is  27.47  minutes  and  the  hydrograph  volume  is  1.73  inches.  When  interpolated  to 
the  specified  computation  interval  (5  minutes,  NMIN  on  the  IT  record),  these  values  become  4,030.04  cfs,  40.00 
minutes  and  1.73  inches.  This  corresponds  to  an  overall  interpolation  error  of  0.7  percent,  which  is  acceptable 
for  most  purposes.  However,  if  the  interpolation  error  becomes  too  large,  it  may  be  reduced  by  reducing  the 
computation  interval,  NMIN. 

As  an  example  of  the  way  the  program  can  be  used  to  evaluate  the  effects  of  future  urbanization,  the 
following  changes  were  made  to  the  parameters  describing  subbasin  1.  Two  overland  flow  strips  were  used 
instead  of  one;  the  impervious  overland  flow  element  was  sized  as  loi=50  feet;  N=0.3;  and  this  impen/ious  flow 
strip  represented  twenty  percent  of  the  subbasin  area.  The  pervious  element’s  representative  flow  length  was 
reduced  from  500  to  400  feet,  and  it  represented  80  percent  of  the  basin.  The  other  parameters  were  held 
constant.  The  results  from  this  run  are  shown  in  Figure  A-5.  The  peak  flow  at  the  basin  outlet  increased  from 
4,030  cfs  for  the  initial  basin  condition  to  4,803  cfs  for  the  increased  urbanization  in  subbasin  1.  The  peak  flow 
occurs  five  minutes  earlier  due  to  the  increased  urbanization. 


Parameter  Sensitivity 

Additional  HEC-1  runs  were  made  with  changes  to  some  of  the  watershed  parameters  to  give  an  indication 
of  the  sensitivity  of  the  modeling  process  to  changes  of  various  magnitudes.  Parameters  were  only  adjusted  for 
subbasin  2.  The  flow  output,  however,  includes  routed  flows  from  the  unmodified  subbasin  1 . 
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COMPUTED  FLOW,  WITH  URBANIZATION 
_  COMPUTED  FLOW,  WITHOUT  URBANIZATION 


Figure  A-5 

Computed  Hydrographs  at  Outlet  with,  and  without  Urbanization 


Analysis  of  Q  Versus  Surface  Roughness 


The  effects  of  variation  of  the  surface  roughness  parameter,  N,  for  overland  flow  is  shown  in 
Table  A-2.  Low  values  of  N  produce  high  flows  and  low  times  for  the  hydrograph  peak.  High  N 
values  produce  lower  peak  flows  and  retard  the  peak.  Note  that  the  high  values  of  N  are  outside  the 
acceptable  range  shown  in  Figure  3  (page  6),  but  are  used  here  for  trend  extrapolation  only.  An 
order  of  magnitude  change  in  N  values  results  in  a  27  percent  change  in  the  flow  rate.  This  variation 
of  results  with  changes  in  N  is  shown  graphically  in  Figure  A-6. 


Table  A-2 

Sensitivity  to  Changes  In  Surface  Roughness 


Case 

Subbasin  2 

Ni 

Percent  Original 

Ooeak  Total  (Cfs) 

TP  (min) 

1 

.075 

.10 

25 

4218 

25 

2 

.15 

.20 

50 

4173 

30 

3 

.23 

.30 

75 

4110 

35 

4 

.30 

.40 

100 

4030 

40 

5 

.45 

.60 

150 

3826 

40 

6 

.60 

.80 

200 

3534 

45 

7 

.75 

1.0 

250 

3308 

45 

Figure  A-6 

Graphical  Representation  of  Sensitivity  to  Changes  In  Surface  Roughness 
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Analysis  of  Q  Versus  Overland  Flow  Length 

As  shown  in  Table  A-3,  increasing  the  overland  flow  plane  reduces  the  computed  peak  flow  and 
delays  the  hydrograph  peak.  Using  short  lengths  has  the  opposite  effect.  The  results  are  shown 
graphically  in  Figure  A-7. 


Table  A-3 

Sensitivity  to  Changes  In  Overland  Flow  Length 


Case 

^trlp  2  In  SubbasIn  2 

®peak  Totai 
(cfs) 

TP  (min) 

1 

80 

4248 

30 

2 

120 

4129 

35 

3 

160 

4068 

35 

4 

180 

4030 

40 

5 

220 

3996 

40 

6 

260 

3858 

40 

7 

300 

3696 

45 

1_ 

oven  I  anci  T  I  ow  i  engt  hi 


Figure  A-7 

Graphicai  Representation  of  Sensitivity  to  Changes  in  Overiand  Fiow  Length 
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Analysis  of  Q  Versus  NDXMIN 

As  mentioned  in  Chapter  1,  NDXMIN  is  an  optional  parameter  used  to  force  the  HEC-1  program  to 
use  a  smaller  distance  step.  The  distance  step  used  in  the  analysis,  Ax,  is  defined  as  the  minimum  of 
the  calculated  Ax  and  L/NDXMIN.  This  NDXMIN  parameter  was  varied  between  acceptable  bounds  (5 
<  NDXMIN  <  50  for  overland  flow  planes)  for  the  overland  flow  plane  in  subbasin  1  to  determine  its 
effect  on  the  resulting  peak  flow  rate.  The  results  are  shown  in  Table  A-4. 

The  first  case  in  the  table  is  from  an  HEC-1  run  with  NDXMIN  unspecified.  Note  how  the  peak 
flow  rate  increases  with  NDXMIN  when  NDXMIN  is  less  than  20.  However,  any  increase  of  NDXMIN 
above  this  amount  does  not  appreciably  increase  the  peak  flow  rate.  Thus,  the  solution  with  NDXMIN 
equal  to  20  is  a  good  compromise  between  accuracy  and  computational  efficiency.  Although  some 
accuracy  was  gained  by  increasing  NDXMIN,  the  interpolation  errors  reported  in  Table  A-4  are 
relatively  low.  Thus,  the  default  Ax  selection  would  most  likely  be  adequate. 


Table  A-4 

Sensitivity  to  Changes  in  NDXMIN  for  Subbasin  1 


Case 

OXMIN 

Q,  (cfs) 

%E, 

1 

5 

2138 

1.9 

2 

10 

2170 

1.1 

3 

15 

2190 

0.8 

4 

20 

2199 

0.7 

5 

25 

2205 

0.6 

6 

50 

2206 

0.6 

Another  kinematic  wave  routing  parameter  which  can  be  adjusted  is  the  slope.  However,  both  the 
slope,  S^,  and  the  roughness,  N,  are  combined  in  the  parameter  a^,  and  either  may  be  adjusted  to 
produce  the  desired  effect.  However,  as  indicated  by  Equation  4  (page  34),  is  a  function  of 
therefore,  changing  has  less  relative  effect  than  changing  the  roughness.  For  example,  to  produce 

the  same  change  in  as  caused  by  reducing  N  by  a  factor  of  two,  the  slope,  S^,  would  have  to  be 

four  times  as  large. 
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Appendix  B 

Sample  Data  Tabulation  Form 


Subbasin  1 

1 .  Overland  Flow  Strip  1 

2.  Overland  Flow  Strip  2 

3.  Collector  Channel 

4.  Main  Channel 


Subbasin  2 

1 .  Overland  Flow  Strip  1 

2.  Overland  Flow  Strip  2 

3.  Collector  Channel 

4.  Main  Channel 


Subbasin  3 

1 .  Overland  Flow  Strip  1 

2.  Overland  Flow  Strip  2 

3.  Collector  Channel 

4.  Main  Channel 


Subbasin  4 

1 .  Overland  Flow  Strip  1 

2.  Overland  Flow  Strip  2 

3.  Collector  Channel 

4.  Main  Channel 


Subbasin  5 

1 .  Overland  Flow  Strip  1 

2.  Overland  Flow  Strip  2 

3.  Collector  Channel 

4.  Main  Channel 


Subbasin  6 

1 .  Overland  Flow  Strip  1 

2.  Overland  Flow  Strip  2 

3.  Collector  Channel 

4.  Main  Channel 


Percent 
Total 
Subbasl 
n  Area 


Cftannal  Channel 
Roughness  Shape  Size  (ft) 
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